SUBCRITICAL FLOW PAST SLENDER BODIES BY A 
TRANSONIC INTEGRAL EQUATION METHOD 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

J. P. AGARWAL 


to the 

DEPARTMENT OF AERONAUTICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

SEPTEMBER, 1979 



OIRIii'IGA'iE 


f2his is to certify that the work » ' SUBORIlDiGAIi 
PASS SliMPER BOIJIES BY A IRANSOKC IHIBGRAL BQUAICION lAETHOB* » 
has "been cjarried out under my supervision and has not been 
submitted eisewhere fca; a degree. 


SeptemHar, 19’?9. 


Ni^ 

(N. L.'^AROKui) 

Professor 

Beparimient of Aeronautical Engineering 
Indian Institute of technology, Kanpur » 


I-OSrGRAIKMTE ■ 

Th,s ,h«« l.„ 

oi iccniioio-v (Al 1'.-^ h ) 

, ‘'-I ^'-cordance vv;di (iic 
I 'e^-Haiious af ,|,e Indian 
I fnsevuic o/ rcchnol, 


Djird. 





ACKNOWLEDGE M ENTS 


I express my sincere thanks to Dr. N. L. Arora 
for his guidance and encouragement. 

I thank Mr. S.S, Pethkar for his accurate typing 
and Mr. Tiwari for cyclostyling. 

Finally, I wish to express my indebtedness to ray 
friends who helped me a lot at different stages of the 
work. 


September 1979. 


J. P. AGAMAL 



SABLE OE GONSEISISS 


Page 

iii 


ABSmOS 

LIS I OE SYMBOLS 

LISS OE EIGFURES 

GMSER 1 UrSROnJOSIOB 

1.1 General 

1.2 Literature Survey 

1.3 Present work 

GHiiPSER 2 PEOBLB^ SPEGIEICASION 

2.1 Introduction 

2 .2 Axisymetric flow 

2.3 Swo- dimensional flow. 


iv 

vi 


1 

3 

7 


8 

8 

12 


OHAPSER 3 INSEGRilL EQUATION MESHDI) OE SOLUTION 

3 .1 Introduction 

3.2 Axisyicmetric flow 

3 ,2 .1 Volume integral 

3.2.2 Surface integral 

3,3 Two-dimensional flow 
5 ,3 ,1 Surface integral 

3.3.2 Line integral 


14 

14 

17 

18 
23 
26 
27 


CHAPTER 4 NIBIEEIGAL SGHEMB 


4 .1 Introduction 

4.2 Axisymmetric flow 

4.3 Two-dimensional flow 


31 

31 

36 


OHAPIER 5 RESULTS AND HSGUSSION 


5.1 Introduction 

5.2 ilxi symmetric flow 

5.3 Twcwdimensional flow 

5.4 Gonolusion. 


40 

40 

43 

44 



ii 


REE'ERMOBS 

45 

iiPRMIEX 

A 

48 

i^PENECX 

B 

51 

APPERJT.X 

0 

53 

itt>PENDIX 

D 

55 

iPPENEIX 

E ■ 

58 

APPENDIX 

E 




ABSmOg 


She sutocritioal transonic flow past axisymotric and 
two-diaensional slender "bouies has "been studied in the 
present work, ghe integral equation oppa^oaoh has heen 
adopted. For this nonlinear transonic small disturhanoe 
equation is transformed into an integral equation using 
Green’s theorem. She solution fcr the axial perturbation 
velocity at a general point in the flow field is sought in 
the meridian plane lay dividing the region of integration 
into reotangular elements, wherein the axial perturbation 
velocity is taken uniform. Ihus an iterative scheme is 
developed of whioh convergence is sought. She results are 
oomputed for a parabolio arc of revolution, parabolic 
airfoil and HADiWOOlS airfoil. She results obtained are 
coii|)ared with experimental results and/or results obtained 
by advanced relaxation methods. 
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FR 

(Vo) 

1 

M 
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x,r,e 


3t,r,e 


x,y,z 

•M 

x,y,z 

Y 

? . P /V 



0 


Pressure coeffioient 

Gcmplete elliptic integral of second kind 
Ocmplete elliptic integral of second kind 
Fineness ratio of paralolio arc of revolution 
tunael half height to body ; chord ratio* 
Laplacian operator 
Mach nuciher 
Porosity parameter 
normalised cylindrical co-ordinates 
^transformed cylindrical co-ordinates 
normalised cartesian co-ordinates 
lEr ansf ormed cartesian co-ordinates 
Ratio of specific heats ( = 1,1 for air) 
Running co-ordinates in x,?,© directions. 
Running coordinates in 5 , y directions. 
Perturbation velocity in x-direction, 
linear pe..'turbation velocity 
tChickness ratio 

Perturbation velocity potential 
transformed perturbation velocity potential 
Parameter satisfying lap lace equation 
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V 


Xj JC j0 

x,S,9 
XjJf z 
Sfy,z 
5,P,V 

?,c 


Values of variable at different locations 
Partial derivative of variable v/»r.t. x,r,e 
Partial derivative of variable wl^*t, Xj^jS" 
partial derivative of variable w.r.t* x,y, z 
Partial derivative of varialiLe wa:.t. x,y,z 
Partial derivative of variable v/,r*t. 5,p,v 
Partial derivative of variable Vir.r.t.5,? 
Pree stream condition. 
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1« Slender body in axisyrametric flow. 

2 , Tbin airfoil in two-dimensional flow. 

3 . location of grid point in reotangular elements . 
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4. Location of grid points in rectangular elements. 
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5. IbLstribution of Gp on parabolic arc of revolution 
(M^ =* 0 .9, FR - lO, Free air case) . 

6. Distribution of Gp for parabolic arc of revolution. 

(M^ =0 .9, RR = lO, Rree air case) . 

7. Distribution of Op for parabolic arc of revolution. 

(M^ ss 0.9, RR es lo) wall interference case with 
porosity parameter 0.77 and wall radius 1.17. 

8. Distribution of Op on parabolic arc of revolution 
(M^ s= 0 .85, RR =s lo, wall radius 1.17) . 

9 . Distribution of on parabolic arc of revolution 

(M^ = 0.9, RR = lo, Solid wall at 1.17) . 

10 . Distribution of Op on parabolic arc of revolution 

.9, RR « lo, waJ-l -- at 1.17). 

11. Distribution of on parabolic arc of revolution in 

free air at supercritical shock free mach numbers (RR = lo) . 

12. Distribution of on gonoral parabolic- body of revolution 
in free air (M^ - 0.9, RR « lo) . 
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13 . Distribution of Op on HAGA-OOlS airfoil at zero 
incidence in free air (M^ » 0 *' 72 ) . 

14, Distribution of Op on 6 percent thick parabolic 
airfoil at zero incidence (M^ £= 0 .82c ) 
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airfoil at zero incidence with wall interference. 

(M^ = 0 »825, h/c =5 2 .0) . 

16. Distribution of Op on 6 percent thick parabolic 

airfoil at zero incidence at Mach number 0 ♦7 and 0*8 
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Despite imoh advances in aerodynamic research for 
flight regime lying far beyond the so called sonic barrier^ 
there still remains some problems connected with the 
phenomenon that occur in the flight range straddling the 
speed of sound. Consequently this range » transonic flow 
range •• continues to engage substantial research effort 
and gained impetus due to possible development of minimum 
drag transonic cruise configuration for transport aircraft 

To simplify the study of transonic flow and for ease 
of building the mathematical model, following assumptions 
are made t 

i) the fluid is ideal and perfect gas 5 
il) the medium is continuum and homogeneous 5 
lii) the flow is steady and irrotational 
and iv) the body considered is slender. 

Under these conditions with perturbation velocity 
(0 i 0 i ) being quite small the equation for perturbation 
velocity potential is given by Prandtl-Kjlauert equation s 


(1 » ) 0^ **' 0yy ^2Z 


0 
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The above equation governs subsonic and supersonic 

flow quite reasonably# However since in the transonic 

range I'loo'* ^ "the coefficient of h.H.S. becomes very 

small hence an order of magnitude analysis of gas dynamics 

equation necessitates retention of some non-linear terms 

27 

also. The final equation takes the form t 

(1 - Moo) ^ 2 Z “ OO ( Y+1)^25- ^ 205 * 

The above equation in cylindrical co-ordinates is : 

(1- ^rr *"2 ^00 ~ ^ Y+l) 0^ 0y^ 

r 

This small perturbation equation is equally valid 
for transonic, subsonic and supersonic flows. In the x-r 
plane, it is a mixed hyperbolic parabolic and elliptic type 
nonlinear equation i,e, for : 


Ql - Ho? (y+1);0X'^ 


> 0 elliptic (subcritical flow) 

= 0 parabolic (critical flow) 

< 0 hyperbolic (supercritical flow) 


The shock free supercritical case is of particular interest 
as it gives minimum drag for the body. 


To solve the transonic flow equation many methods 
have been tried of which relaxation method is quite standard 
method but needs a lot of labour. An equally good technique 
is integral equation approach. It is less time consuming 
and attempts are being made to make it of similar accuracy 
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as tjiat of relEucatictti method* A brief survey of integral 
equation method is presented in the following section. 

1»2 Idteratune Sur^y i 

Integral equation approach is oie of the earliest 
successful method for solving the direct transonic flow 
problem. !Dhe main idea is to convert the transonic flow 
equation into an Integral equation using Green’s theorem. 

The solution to the integral equation is obtained by 
iterative technique after certain modifications have been 
incorporated into the equation so obtained. 

Among the most noteworthy earliest atteiai^ts for the 

solution of transonic flow problem by integral equation 

method are three related studies that deserve special 

mention. These include those carried out by Oswatitsch, 

Gullstrand and Spreiter and Alksne, and hence the approach 

1 

adopted has been termed as OGSA approach by Ferrari « 

2 

Oswatitsch deduced a ncn-linear integral equation over the 
whole plane flow at zero lift and made approximate calcu~ 
laticaas for some aerofoil in the lower transonio range. The 
work was later extended and iHustrated with applications by 
Gullstr^d ■ • more recent exposition of the method is 

Q 

given by Spreiter and Alksne , who worked out the method 
more thoroughly* Their results are fairly good for the 
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subcriijioal flow around circular arc airfoil at zero 
incid^oe * Heaslet and Spreiter^ applied the integral 
eq,uati< 3 iQ. njethod to transonic flow around slender wing and 
bodies of revolution* Special, attention was given to 
conditions resulting from shockwaves. Results are obtained 
for cme cylinders, wings and wi3ag body combinations at a 
free stream Mach no**«l and compared with experimental 
results • 

N^rstrad"^” extended the OGSA integral equation 

method to lifting and non-liftirg transonic flow problems, 

including the flow with shock discontinuities. She integral 

equation is transformed into a set of non-linear algebraic 

equations* (Do ensure a unique solution for supercritical 

shock free flows the method of parametric differentiation 

has been applied. His results indicate generally good 

agreement with experimental results and otherwise- accurate 
those 

results IjOce^/obtained by finite difference techniques. 

(Bhe integral equation method to solve the problem of 
high subsonic flow past a steady two-dimensional airfoil has 
also be^ used by Nixon and Hanoock^^, (Dhey have confuted 
the shook free two-dimensional flows around the lifting and 
non-lifting airfoils by solving the integral equations 
approximately, (Dhe results for two test cases show close 
agreen^nt with more exact numerical results of Sells^, 
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The method has also been extended hy Nixon to an airfoil 
oscillating at. low frequency in high subsonic flow. 

A feature common to all above integral equation methods 
is that the field integral is not evaluated accurately. In 
these methods using a functional relationship for velocity 
distribution, the double (surface) integral is reduced to 
a line integral in the plane of the airfoil. In the extended 
integral equation method of Nixon an alternative means of 
evaluating the field integral is developed. Ihe flow field 
is divided into a number of streamwise strips and the trans- 
verse variation of perturbation velocities across each of - 
these strips is approximated by an interpolation function 
in terms of values on the strip edges, Ihe field integral 
is then reduced to a line integral \'rhich is in turn evaluated '*• 
by quadrature. Ihe results obtained by using this method 
have been compared with those of the standard integral 
equation method (See Nixon ), It is shown that this method 
gives considerably improved results over the earlier integral 
equation methods. 

Recently Ogana^^ obtained results for two-dimensional 
parabolic arc airfoil evaluating the field integral quite 
accurately. He divided the region of integration into 
closed rectangular elements, and considered velocity to be 
constant in eaoh element 5 the integral equation is then reduced 



6 


to matrix equation solvable by suitable iterative scheme for 
subcritical flows • 

fhe influence of wall interference effects on models 

tested in transonic wind timnels has been an open question 

for many years* fhere has been an increasitg need for the 

largest model possible. She relaxation method has been 

applied to the problem of wall interference effects. Bailey^® 

treated the slender bodies of revoluticaa at transonic 

speeds in the presence of wind-tunnel wall boundaries, while 

Murman^^ and Kaopraynski^^ have considered interference 

effects on two-dimensional airfoil. Newman and Klucsker^^ 

have extended relaxation method to determiiB interference 

effects of transonic flow over finite lifting wings in 

rectangular tunnels. Murman, Bailey ani Johnson^ have 

given an upto date conputer program to solve the transonic 

small disturbance equation for two dimensional flow for 

liftiig airfoil which takes free air case as well as various 

wind tunnel wall conditions. Ihey used finite difference 

formulation using iterative successive line over relaxation 

(SIOR) algorithm. An integral equation method has also 

25 

been used by Kraff for wind tunnel boundary interference 
at transonic speeds on thin airfoil in two-dimensional wini 
tunnels. However, field integral is reduced to a line 
integral over the airfoil hy the use of arbitrary approximate 
function and the results are obtained on an airfoil. 



The integral equation method applied to two- 
dimensional flow has proved to givd good results. However, 
the attempts to solve the problem of mixed transonic flow 
past bodies of revolution by integral equation method have 
been dismal. 

In present work an integral equation suitable for 
axisymraetric flovr past slender bodies of revolution at 
transonic speeds is obtained in the free air and in the 
presence of porous wall boundaries. The numerical solution 
to integral equation is sought in the meridian plane by- 
dividing the region of integration into rectangular 
elements wherein the velocity can be considered uniform. 

The pressure field is computed on and away from the slender 
body including force and aft body effects. An application 
to parabolic bodies of revolution with and without sting is 
presented at subsonic Mach numbers yielding suhcritical 
and slightly supercritical shock free flows. 

Similar analysis has been carx’ied out for two- 
dimensional flows also . Parabolic airfoil has been 
considered in the free air and in the presence of porous 
wall. An attempt is also made to determine pi’essure 
field for NAOA 0012 airfoil at zero incidence. 
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PROBM aPBQI F IQAglOK 
2,*X Introduction ; 

Most of the aerodynamic problems are related with 
finding pressure field over two-dimensional or axisymmetric 
bodies, so that one can predict coefficients of lift, drag 
or mcment over a body in certain situations. Efere we are 
interested in finding pressure field on and away from body, 
fhe axisymmetrio problem is considered first and then the 
flow past two dinensional non lifting airfoil is treated. 

2 ,2 Axisynmetric Plow s 

Consider the inviscid, compressible flow about a 
slender, pointed body of revolution in cylindrical cos® 
ordinates x and r (made dimensionless with respect to the 
length of the body) with x axis parallel to the body axis 
and free stream velocity (see fig, l) , Assume shook 
free flow so there is no vorticity involved. Under this 
assumption a velocity potential can be used to calculate 
the flow field. In particular a perturbation potential 
can be defined such that the perturbation velocities, made 
dimensionless with respect to the free stream velocity, are 
h- — 02 J- parallel to the x~>axis and v = 0^ parallel to r— axis. 
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The governing differential equation for transonic flow 

given by slender body theory can be written as^*^ 

2 

where is free stream mach number and y is the ratio 

of specific heats. 

The term associated with 0^ mates transonic equation 
both nonlinear and of mixed elliptic-hyperbolic type. The 
mixed character of the flow field here may occur with 
local supersonic regions embedded in a subsonic flow. 

To coD^lete the specification of problem* boundary 
conditions must be given at the body and in the outer flow. 
The flow tangency condition at the body surface given by 
the ficst-order slender body approximation is that near 
the body axis^"^; 


lim- r 0 « R 
r o dx 



( 2 . 2 ) 


where r = R(x) defines body 

S(x) defines cross sectional area of the body and 


prime denotes derivative of function w.r.t. its 
argument • 

For a body in free air the perturbation velocities 
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vanish at ini’ ini tyj 

02; 0 ^ ( 2 *3 ) 

which is satisfied by setting 0 equals a constant, 
say zero, at infinity. 

In addition to the free air boundary condition wall 
boundary condition may be given to approximate inviscid 
flow in porous wall wind tunnel test sections, Ihese 
boundary conditions can be used to illustrate wall-induced 
interference effects. Although the formulation is strictly 
valid for a circular test section, the results may be 
compared with a square test section of equal cross sectional 
area because the effects of wall Interference at tie centre 
of tunnel should be relatively insensitive to the actual 
wall shape • 

Ihe average boundary condition for a porous v/all, as 

Pfi 

derived by Goodman ^ follows from Darcy »s law for slow 
viscous flow through a porous medium. It is assumed that 
the average velocity normal to the wall is proportional to 
the pressure difference across the wall, which is a linearized 
approximation of Darcy* s law for a thin wall, and pressure 
outside the wall is free stream pressure . With the wall 
parallel to x-axis the porous wall boundary condition becomes: 

tp02^ ~ 0 at the v/all r = rw (2,4) 
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The quanti1;y p is a porosity parameter or Reynolds Number 


of the porous medium, defined by t 
P 


U P k 

oo oo 


w t 


(2.5) 


where 

free stream velocity 
free stream density 
k permeability of porous medium 

y viscosity of air 

t tunnel wall thickness 


k is determined by the structure of the porous 

26 

medium and must be found experimentally . 


As porosity parameter vanishes i.e. p = 0 equation 

(2.4) gives boundary condition for solid v/all and the 
porosity parameter becomes very large' i.e, p equation 

(2.4) approaches the open jet boundary. Although truly 
porous walls are seldom used for transonic wind tunnels, 
the porous wall approximation is useful mean boundary 
condition for perforated walls. 

Finally, the slender body approximation for the 

p Jn? 

pressure coefficient at points near the body is given by 

= (g.6) 

and cn the body surface 

Op = - -(dE/te)^ 


( 2 .'!’) 
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1«3 (Pwo-Dimensicmal How j 

Here Cartesian coordinates (made dimensionless 
w .r .t . the length of airfoil) has he on chosen to describe 
the in viscid, compressible, shock free flo?/ over a thin 
symmetric airfoil, 3}he x-axis has been taken to ie parallel 
to airfoil plane and the free stream velocity (see fig,2), 
Ihe perturbation velocities made dimensionless w.r.t, the 
free stream velocity U , are u = parallel to the x-axis 
and V a parallel to the y-axis . 

IBhe governing equation for transonic flow 
given by thin airfoil theory can be written as^*^ 

(1- + 0yy = 0^ +1) 


where symbols have same meaning as in pre-vious sectiont 
She boundary condition for thin airfoil is 

■u 27 

given by 

lim- 0^(x, +0) a (2,9) 

y *♦ + 0 dx 


where y = 


' S+(x) defines upper airfoil profile 
Y-(x) defines lower airfoil profile 

S'or symmetric airfoil Y+(x) = -Y— (x) • 


For the airfoil in free air perturbation 
■velocity vanish at infinity which is satisfied by setting 


*^x » ^y 


y^ 


( 2 . 10 ) 


*-► OQ 



For wall interference case the effect 


porous wall, parallel to free stream can le taken iDy 

?6 

applying porous wall boundary condition given by 
0^ + p 0^ = 0 at y = + y^Y 

where + sign is for upper wall, - sign is for lower wall, 
and p is defined by (2.5), 

I he thin airfoil approximation for pressure 
coefficient for the point near the airfoil is given by 

Cp =* 2/ynliii >2 _ ^2jj 

M X y 

( 2 . 12 ) 



OHAPaiBR 3 


IHiEEGRAIi EQUAIION lilSgHQD OF SOLUTION 

5 ,1 introdUGliiQn t 

In this chapter the governing nonlinear 
partial differential equation of transonic flow is trans- 
formed into an integral equation using Green’s theorem for 
subcritioal and shook free supercritical speeds (M^ < 1) . 
Axisymmetr io and non-lifting two dimensional cases have 
been dealt with* !Che axisymnetric case will be treated 
first • 

3 «2 Axisymmetrjc flow : 

So apply the Green’s Iheorom to the axisymnetric 
transonic flow equation (2,1) one stretohis the coordinate 
system suitably. Here we make the following transformations: 

x=x; 0 = 

Tf ' 2 

where p 1 - 

K = D 4 

Now the goveming transonic flow equation for axisymnetric 
flow (eq« 2,1) can be expressed as 
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^ XX 


(3.2) 


Similarly boundary conditions equations (2,2), (2,3) and 
(2.4) beccme 

0 at infinity (3 ,4) 

and for wall interference case at the wall r =s r^^- s 

h + = 0 (3 . 5 ) 

Now we shall apply Green’s theorem to the equation (3,2), 

The theorem states that the following relation holds between 
any two arbitrary functions f and g having continuous first 
and second derivative in the volume V bounded by the surface 

S s- 


fffZsHt) ~ fL(g)nw = //(^ H - s H ) as 


V 


8n 


(3 ,6) 


where n is inward normal drawn to the surface S and 1 is 
Laplacian operator in cylindrical coordinates i.e, j- 


li s 


a ^ i ^ + a 

ar p 9p 9p 


2 + “2 


1 

( 

p 


a v" 


here (^ are running coordinates of domain ? and S, 
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!Bo obtain the solution for 0 one chosess 


f (C rp»v) - ^(5»p»v) 

= 0 (Op ) because of axisynmetric 

nature of 0 

sC? »p>v) = 4) (c»p>v) 


where ifi satisfies L( tfi ) = o which gives the fundamental 
solutions 



where H « [;(x «C )^+ r% 


(3.7) 


2 rp 


cos(0-.v ) 


Vg 

a 


Shen Eq.uation (3 .6) gives the integral equation s 

///♦iC?) dv = // (^ II - ^ M )as (3 .8) 

V s 

mm 

Since rj/ and 0 do not have continuous first 
and second derivatives everywhere in the flow field hence 
we shall investigate to what domain above relation holds, 
fhe volume is (unhatohed portion in fig. 1) infinite volume 
Vg, (enclosed by wall if wall is present) excluding (i) a 
spherical point cavity at (c = x ,p = r , v = 0 ) where 
is singular and (ii), the volume which the body 

covers, fhs surface S bounding the volume V consists of 
i) surface Sg, at «> for free air case or at wall 
if wall is present ; 
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ii) surface Sp of spherical point cavity at 
( 5 = X ,, p = r , V = 9 ) 
and iii) surface Sg hounded by the body. 

Now we consider the volume a-id surface 
integral separately, 

3 ,2 .1 Yolume Integral s 

Prcm eq.uation (3,2) 

1 m 

where u » ^ 

hence the volume integral in e equation (3,8) becomes! 

iy= I SSI f l^av (3.9) 

here T = Vjj, 

Ihe cavity volume Vp can be included in Y 
as it gives negligible contribution in the limit to (3 ,9) 
as shown below, 

let us assume the radius of spherical cavity 
to be G I e “• 0 , we get 

4 3 

. dV « Yp = y tt e 
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.^2 


= finite ss G^(say) 


r 1 

J7 


Iff 


r-2 


au 


Vt 




dV 


* I 

= |0l®^ 


» 0 


as G 


0 


hence we get 

V = Vj - Vfi 

Now integrating equation (3 .9) ty parts in ^ direction 
with boundary condition (3 .4) of perturbation •velocity 
vanishing at ^ = + “* one gets 

ly = - I 

5 .2 Surface Integral s 

i) Integral over i 

For free air case the integral over the 

oylinarioal surface S„ at infinity can be e:! 5 >bessea as 



«> Sin „ ^ 

w f “ 'I' ^ p dv d^ 

lim p » 
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This inliegralf with 41 given, hy equaijion (3»7) and wiiJh 
the assumption that ^ G >■ 0 at infinity, 


vanishes • 


ffor wall interference case 'uhe integral over 


the finite cylindrical wall surface can he expressed as 


Iq = Iq *s , 
S, 


•• Sf (^4> — ^ ^w 

S.„ P ^ 


where r^ is the radius of cylindrical v/all 

«• 2ic _ 

or le = 5„ / / (♦?„ - 0'1'J „ ■1'’ 45 


- d w ^ 

^ ^ o 


P -r 
P “^W 


j^pplying porous wall boundary ocsndition (3 #5) 


we get 


«> Sic 




P “^w 


After simplification eq.uation (3 .11) leads to 
following line integral (see Appendix A-I) 


Iq ^ J 


/"t:(A+ B) i (C , - 0 p- . 5w)ia5(S.i2) 


where 




r^ - ^ E(it/2 y 

% (a-h) -fa 


1 P(ti;/2, k) - B(V2, k) 

Si Ya 



20 


0 = Ee 

a =: )^ + (^ + ■>3 = 4 r r^, 

k = f B7a 


(3,13) 


E and i* being the complete elliptic integral of second 
and first kind respectively. 

Hence surface integral over can be expressed as 



where 


« 6 

6 


/ [[[.(A+B) 5 ) r^) ■* G o ^ ( 5 , r^) 3^5 

,oo r ^ 

(3 ,14) 


0 for free air case 

1 for wall interference case 


ii) Integral over Spj 

She point cavity is the sphere of radius G» 
G «> around the point ( 5 = x , p = r, v = e) the 

integral over Sp can be written as 



= C I - ^ S 

« illa|4 « -A— M 3 4m: 

^ 4 icR^ W aR 

s S(x, r) - G ^ (since here I^g) 

eR 

sa -* .? (x, r) G ** 0 


(3,15) 




For suf ficiently slender body we approximates 


lim 
P -*o 


d, ^ 
an ■" ap 


hence we gelis 




p <3 p d 5 


(3,16) 


where Z is the non dimensional length of the body • 

— — 27 

In the limit p o, 0 is of the form 


0 « S*(5 ) P + g( C ) 

Also from (3 .7) j 



0 


% 


r oos(e~v ) 

- T”)^ 4- 


5^/2 


(3 .17) 


Hence ih the limit as P**0, Pj2^’l'p=0 and the integral 
(3 ,16) reduces to 



} ^27t 

-S f ^ 

0 0 . 

lim p -0 


p do d C 

ap 


Applying the tangency boundary condition (3.3)* 

I , » /f ^ s:ki^ dv dc 

0-0 Z% 

lim P o 
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^ — 


’ (x~^ ) + r 


(3 .18) 


Now using (3,3X)), (3.14), (3.15) and (3.18) in Green»s 
Iheorem (3,8) we finally obtain, a non-l:inear integral 
eq.uation for 0 

^(x, r) ^ ^ ■) dr 

t; -tr*" ^ 




+ r 


^ jL ^ (e » ^w) - ^ ^ 

(3.19) 


+ I If I (c ,p )dv 


p ^ v» 

V 5 c 

Now to obtain an equation for the axial perturbation 

velocity we differentiate equation (3 ,2o) w.r.t. x 
(see Appendix A.-1I) and we get $ 

Z 


u(x, r) 5s 


o | 2 (x -e )^+ 52^372 


d? 


VA# 

+ 8/[:A+B + D:]u(e, rjdj . 


•.^00 

I 5S 


, P X - € 

where N = r ^ • ' ■ a 


(3 .20) 




and A, B' stand for same notations as defined in (3,13). 
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!Ehe equation (3*E0) can be rewritten as 

•* M oo 

u(5,r) = Uj.(S, r) + (A+BfD) u(C ,r^) d? 

1 

2 '^ 1 / X ^ > p ) 


V 


(3 .21) 


where Uj^(x, r) « ^2 Uj^(x, r) 

P 

n - 1 r s‘(e ) (x-e ) . r 

Uj^ being the linear value of perturbation velocity 


3 .5 !Bwo~ Dimensional Plow 

In previous section we derived axial velocity 
perturbation for axissnumetric flow. In the same manner 
the axial velocity perturbation in two dimensional flow 
over thin airfoil at zero incidence has been derived in 
this section. 


We first introduce the similar transf ormaticn 


• 0 • 


xaxj y~py?^ 



(3.22) 


where p s= "^1 « K = (Y + 1) Klf, i 

GO 4MW. 

Using (3,22) in the governing transonic flow equation (2.8) 
and boundary conditions (2.9), (2.1o) and (2.11), we obtain 




(3 .23) 
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(55, oi) 


aY^(x) 

dx 


where Y (x) »= — p Y^(x) 

t- R + 


-* 0 at infinity 


(3^4) 


(3.25) 


and for the wall interference case, at the wall y = + 

_ p- 

% = + ,p 0^^ (3.26) 

For two-dimensional case the Green’s theorem which relates 
the surface integral to line integral is : 

// dS = / (f II - g II ) dO (3.27) 

S G 

where f and g are arbitrary function having continuous 
first and second derivative in surface S = S( ? , ? ) 
bounded by curve 0 = 0( C » C ) and n is the inward drawn 
normal to the boundary • L is the laplacian operator given 
by 

L ^^4. 9^,. 

3 5 "^ d C ^ 

So investigate the solution for 0 , we ohooset 
f W ^ ( 5 , C ) 


and g « ,|) 
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where \p satisfies the Laplace equation L(i|; ) = o. She 
fundamental solution of this equation is : 

t|; « On R (3,28) 

where R = 

Hence equation (3,27) together with equation (3,23) gives us: 
//♦ I = If) dO (3.E9) 

Hov\r we shall investigate the surface S and hounding curveO 
v/here 'P and 0 have continuous first and second derivative 
keeping in view the boundary ccaaditi(m. 

Such surface is infinite plane S^, (enclosed by wall 
if wall is present) excluding (i) a circular point cavity 
Sp at(5=x,C ®y) (ii) the airfoil surface Sp . 

(see fig , 2) , 

$he cuT’ve 0 consists of (i) curve around infinite 
plane (or wall if wall is present) (ii) Qurve Op around 
ciroular cavity and (iii) curve Op around airfoil surface . 
(see fig , 2) • 

How we shall analyse the surface and line integrals 
of equation (3 ,29) . 
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5 «3 Surface Integral; 


given iiy 


Frcm equaticai (3 ,29) integral over surface is 

Is = //'" 4S 

o 

= i " 0 ^ ^ ^ *30 ) 


}aere S = S^, - Sp - 


She cavity surface Sp can be included in S 
as it gives zero contribution in the limit to (3,3o) as 
shov/n below# 

Let radius of circular cavity be e, G ** 0, 


we get-5 

dS = 4u: 

u^ =! finite = Cg(say) 

= I (In e) 


(since here R = G) 


=5 0 as G 0 - 


hence we gets 

S = Sg, - Sp 

How integrating ec^uation (3,30) by parts in 5 direction 
with the boundary condition of perturbation velocity 
vanishing at infinity we get s 

I, = - I ff ♦5 ^ as 
s 2 g 


(3,31) 
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3,3^ 


i) Integral over 0^ ; 

Since all perturbation vanivsh at infinity henoe 
for free air case integral over Ofji gives zero contribution 
while for wall interference case integral over gives a 
finite value given by (see Appendix B) 

oo 

~ ^ jL ^ 11 ♦ ( c = y^) + 1' ( ? =-yw) U 


- / T'l’c («= yw> -1'!; «} « (S^2) 

0 for free air 

where $ s 

1 for wall interference case 
y at upper wall 

and 5 ss 

— y^ at lower wall 


ii) Integral over : 

Integral over curve Op is given as follows 

\ ^ n - ♦ ^n) iO 

*= y) Ir - I/) ^ I] 2tcG 

dR aR 

where 6, e o is radius of circular cavity Gp . 
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I = 2it ? (x, y) _ 2n G In 6 

^ 8R 

(since here R= G) 

« 2it ^(5c, y) as e - 0 (3,33) 


iii) In-feegral over Gg s 

Integral over hody curve 0^ is given by ; 

B 

“-/cm G5-5)®4. 

« A ^ 2 dC 

(c -x)^ + ( c -y)^ c »o 

= ^ (x, O"^) ~ 0^ (x, 0“) 

= ? (5j o"^) - 0(x, 0") 

= 0 for symmetric airfoil at zero incidence 


here ^ 0 
A ? 


from boundary condition (3 ,24) we get: 




( 5 ) 
*5 



dY_( 5 ) 
dC 

for symmetric airfoil at zero 
incidence , 


which leads to 



2 /in C(5 -5)®+ 


dr+(g ) 

d5 


d? (3 ;3 4) 



OomlDining (3,31), (3.32), (3,33) and (3,34) in equation 
(3,29) we get integral equation for the transformed 
velocity potential as 

0(Xi y) - - / :in [; (C - x) + <3. K 

+ fei ( 5 » yw) { I C (^ = yw) 

^oo P 

- f (c = %) - =- yw) H 

- ^ 0^ <3-5 dC (3,35) 


To obtain integral equation for axial con^jonent of pertur- 
bation velocity we differentiate equation (3,35) w,r,t, x. 


Thus 


tl(x, y) « 


K 


f 


^ X - K 


( O 




<3-C 


2 % 


^ P 

(c . yw) { p Ci-j (c=yw) 


+ 1' (« = - Yw) H -l-r ( ? = Yw) 

- k ^ ^ 

GO 

» Uj^(x, y) - Ij/ u(5 .Yw) { f- lli(ij.(? = jfw) 

*— oo P ^ 


or iX(x, y) 
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where (x, y) = Ut.(x, y) 

/- -X 1 X ~ r dY. ( O 

U^U, y) (3.37) 


here Uj^ is linear value of perturbation velocity for 
symmetric airfoil at zero incidence . 



OHAPi m 4 
TOIERIOAl SCHEME 


4 ,1 I: troductlon i 

In -this chapter a numerical scheme is presented 
which has "been used in calculating integrals involved in 
computation, of perturbation velocity, She integrals have 
been suitably expressed in summation form approximating 
velocity perturbation to be constant in small intervals. 

In this manner a set of algebraic equations are formed 
which are solved by method of iteration using suitable 
ccnvergence criteria. 


4 ,2 ; Ax isymmetric Fl ow s 

We consider equation (3,21) for numerical 
evaluation. Uj^(x, r) for general axisymmetric body can be 
calculated numerically as follov/s . 

Ul, (x, 5) = 'i 5 


which on integration by parts give 


Un 










fi ^ P' 

S (x) In 


4x(2, «x) 
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On the body surface r = R(x) or r = pR equation (4 . 1 ) 
for pointed nosed bodies simplifies to 


u,(s, 0) = r.^ 


Jt - X 


I i _ 

+ s (x) In 


4x( S- - x) 


Z 


0 |x -c I 


(4.8) 


Ihe integral in eq* (4,2) vanishes for 5 = x, Ihe simpson’s 
rule can be used to evaluate this integral quite satis- 
factorily. 


For the evaluation of linear perturbation 
velocity away £:om the body, following procedure has been 
adopted} 


n^ 




(X , ?) = ^ s'(Sp / -I - ' 


i=l 


x^-li 


(x - 


d 5 



where x^ - 1 ^ ~ 0 

\ = ^ 


(4 ^ ) 


Here the body length has been divided into n^ numbers of 
intervals of width 21 ^ » 
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Before ocsnsidering equation (3.22) we attach a subscript 
S to X and. P to r . to siiiplify the analysis further. iBhe 
second expression on RHS of equation (3.21) can then be 
e3q)ressed as t 


oo 


6 / (A + B + D) u ( 5 , r ) 

Vv 


n 


" ^ i^l ^iw 


(4 #4) 


whero 


\wSP 






=___ / (A + B + B) d5 


^i"^i 


Here the wall has been dinrided into n number of steps of 
size 25^ and the yelooity at the raid point of each step 
has been assumed constant over that step . 


[Do simplify further we proceed as in Appendix G 

and finally obtain : 


^iwSP ** 


2ic 



tan' 


-1 


(v-s'p) 
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where 


Xg - + 5^ 


Xg - - 6i 


Xg - \ 


^ h % 


X^+ (r^+^p)^ 


0}o evaluate the volume integral of equation (3.21) 
the region of integration has hoen divided in rectangular 
grid in meridian plane, Ihe computational domain is defined 
by grid network as shown in fig. 3. Ihe perturbation 
velocity at the centre of a grid has been taken constant 
throughout that grid. Hence we get the volume integral 
after sinplifioation as (see Appendix D) 


I fff If, - U ,P )<iV 


n m p 

« 2 2 q. 

i=l 3«1 


^ioSP 


(4,6) 


where u. . - 
10 


u(xj^,_rj) 

r .+ h. 
1^3 3 

S' '> 

^3-13 


8lt 

0 r: . _ , K X 


dv dp 


- G(Z2) B(7t/2, kCXg)) ^(Xg)] 


(4:.>?) 
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where G-(X) = 


V r+(rp+ rp2 


T(X) 


tan“^ - tan*”^ ^3“ “ q-n 


k(X) = 




X^+ (Sp + r.)' 


finally using (4.4) and (4.6), equation (3.gl) can he 
expressed in iterative form as 

(n) n (n-1) n m gC^-i) 


J^) 


n 


(n-l) n m 


« Ut 

+ 6 £ 
i=l 

^iwSP 

u • ^ P 

IW 


« g (s 

ay) 



wl^re 

p:i 

II 

U (Xg, 

Hp) 


and 

II 

dp (xg, 

Hp) 



(4.8) 


lo..' convergence of above iterati-ve equation ^ Ugp 
has been added to both sides where X == - = - 2 ®'SPSP^P* 

dUgp 

so that we obtain 


" L ^iwSP '^iw 


(n-l) n m o(M.-i) 

u. + S E a,,^ r.. 

“ i=i 3=1 ^3® 13 


** ^ ®'SPSP '^SP 


II / C ^ 2agpgp Ugp 2 (4,9) 
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Ihe convergence is achieved provided agpgp 


^ 0.5/Ugp 


^ Swo -dimensional Plow s 

Ihe numerical iterative form of the integral equaticxn. 
(3,36) for two-dimensional case can he obtained in the 
similar way as in i 4 ,2 , 


I'or most of the twou. dimensional airfoils analytical 
results are available for linear value of velocity potential, 
fhese results for parabolic airfoils and ITAOA OOXX airfoils 
are presented in Appendix E-XI . 

Ihe wall integral in equation (3,37) is numerically 
ejipressed as 


- ^ f f- ll'i’j (t = %) +>fj (c =-yp U 

(c = y„) + ipjte = -%)y <3£ 


K 6 


n 

E 

i:=l 


^iwSP 



where = u (x^, y^) 


introducing x 
'^iwSP ^ 


- Xg 


s 


^i'^i 


and y = yp 

^ I (c 


= y^p Cc = -yw^H 


(c = y^) +ip^ (c =~ y^)} de 



3*7 


~ If ^ f = yw^ + n 




X.- 6^ 


x.+6^ 


i— *-1 C 1 3 Cq-^ 5 

+ L^!“i _-^ + tan- J—— ;] } 

^v/"*" %> x.«6. 


yw - yt 


1 1 


1 ,p 


" Si~^2p' C 


in 


( yw"* yp ^ "** ^ ^ yw"^yp ^ ^ 

(V yp)®+ ^ " (vyp>^ + 4 


+ 


-.1 - - ^ 


yw " yp 


- tan 


^ ^1 . a.„„-l ^ 


+ tan' 


yw- yp 


yw+ yp 


1 ^1 

tan--^ ~=^i: ■} 

V yp 


where = Xg - + 5^ 


""S - 


Before representing surface integral of equation (3«36) in 
suitable form for numerical evaluation we use the symmetric 
property of 1i( 5 > ? ) for symmetric flow and limit the 
integration region to upper half plane only, as follows? 


" ^ 


1 


/ f 


y^ (jP 


' Xq ) "^yp) 


-Jrp C (§-Xg)^+ (^-yp)^Il 


-.2 




J. 

4n: 


^ p « p (ym is extent of C) 

^ ^ t rr. ! T, \(^ -L 


j” / r -(c-yp) 


•oo o 

- .2 


C(5-%)''+i:?-yp^‘^ 

Cc— Xg) — Cc yp)^ 

c7e-5s)''+ fe*yp)‘=3^ 


+ 


c[(e,Od5 d5 



38 


2 (say) 

Now the region of integration is di-vided into rectangular 
elements, as shown in Pig. 4, She perturbation velocity 
at the centre of each element is assumed constant within 
that element, Ihis gives s 


n m 


2 


13 


^IjSP 

_2 




a. 


IjSE' 


X4+ >4 .0 o 

1 r ^r(5-Xq) - fe-yp)® 

/ / 


V -*-0 • X. — 6 > ^ 
*^3 ^3 ^ i 




(e-xg)^ - (e-f y-^f 
[I ( 5 -Xq) + (c+ yp) I] 


d5 dc 




^ f 


yj+h. 


Xg « c 


Xg - c 


C|'-Sg)^+ C’S-yp)^ (C-Xg )"+(<?* yp)''J 

% 


-“TI2 


V ^i 

dc 


^3^ -3 - 




^3“ 


PM . 

^i*^i 


X. 


;r;F 


L ^ + (? +C5+yp) 

s, Zp 




dC 




-1 4.._-l V 


15. ^ j *b8xi 

y-D- c -~j ^3 ^3 


X-, 


^ .,-1 ^p“ ^ 

+ tan 

r. 


^3 ‘^0 
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• ^ X, Zg 

^ ~1 ^4 ^ ^3 1 Xp . Y, 

•• tan + tan «~ + tan"* -. tan"** wi H 

■^1 = % - 5^3 + ! ■^s = yp + ^3 - ij 

Ig « «. y^ - = Yp + y^ + 

Sims equation (3 ,56) in iterative form becomes 

^(n) ^ ^ (rul) n m g(^^) 

^S9 ** ^SE * ‘^iwSP ^iw ^-iiSP ^ij 

(4.1o) 

where Ugp (xg, yp) 

«Lgp = “!,(% > yp) 



S^JgLS. 

RiSU IigS MI), D ISCU SSION 
5 * 1 Introduotion t- 

In this chapter results of computations are presented 
and discussed. Results are computed at 2 :ero angle of attaolc 
for parabolic arc of revolution, parabolic airfoil and 
NAGA^OOlB airfoil, 2hey are compared with experimental 
results and/or that obtained by relaxation technique, She 
discussion for axisymmetric flow and two dimensional flow 
are dealt with separately, listing of the computer programme 
is given in i^pendix R, Ihe computations were carried out on 
DECU1090 System. 

5,2 Axisymmetric Plow i~ 

Various results for parabolic arc of revolution of 
fineness ratio lo with sting at 0,854 have been obtained at 
different Mach numbers in free air and in the presence of 
wind tunnol wall. Ihe aft and fore effects are also included, 
Ror free air case integration ?/as carried out in the 
range-O *28 £, x 1,16 and 0 r <1,1, ghe streamwise grid 
size of A X = .04 was chosen. Ihe transverse grid size A.r 
in physical plane was .2 except on the body surface where it 
was .1* Overall 37 x 6 rectangular elements were taken. 

Rig . 5 shows distribution of Op on the body at Mach 
numbac o ,9 . It agrees quite satisfactorily with Bailey and 
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experimental results of laylor and McDevitt^^ except in mid 
region* iBhe results for this case are also presented away 

from the body and conpared ?/ith experimental results^® (see 

fig. 6). lEhe convergence is achieved in 4 iterations'. 

Pig* 7 shov/s a comparison between present result with 

porous wall at radius 1.17 and porosity par^mieter 0.77, with 

29 

the experimental result ^ carried out in similar conditions^O 
Ihey are in close agreement exccept near the sting. Ihe 
integration was carried out in the range 0< x < 0.854 and 
0 1 3x < 1.17 without taking fore and aft effects at Mach 
number 0.9. It took 14 iterations to converge which indicates 
the porous wall interaction v/ith flow perturbations. Ihe 
results are plotted on and away from the body. 

Pig. 8 shows a comparison between three oases i) body 
in free air ii) body in presence of porous wall with porosity 
parameter 0*77 and iii) body in presence of solid wall, at 
Mach number 0.85. She axial integration range was - .28 < x 
< 1*16 and the wall v/as taken at radius 1.17. Par conver- 
gence the number of iterations were 3,9 and 5 respectively. 
2he porous v/ePl case do not have much difference over free air 
case. In the solid wall case the flow shows increased accele- 
ratim in the mid region. 

A comparison betv/een free air case and solid wall case 
ail- Mach no *.9 w'ith fore and aft effects is shown in fig* 9. 



More increased acceleration of the flow is observed here in 
mid region which signifies that with increase of Mach 
number the effect of wall interference increases. IDhe 

convergence was achieved in 4 and lo iterations respectively 
fcr the two cases. Porous wall interference calculations 

did not converge in the limit of 15 iterations. However, 

without for and aft effects it converged in 14 iterations 

as mentioned above. Oorr espondingly solid wall and free 

air calculations were done without fore and aft effects 

(see fig, lo) » Ehey took 8 and 4 iterations respectively, 

2he present method works also for shock-free super« 

critical case at Mach numbers 0.950 and 0.955 for free 

air with aft and fore body effects , Oonvergence v/as 

achieved in 1 and lO iterations respectively, She pressure 

distribution on body surface is presented in fig, 11, fhe 

result of Mach number o .950 shows agreoahle matching with 

2Q 

635) erimental results , 

Pig, 12 shows a sarple calculation for general 
parabolic body of revolution having maximum thickness at 
3o percent and 70 percent of body length, Oonvergence was 
achieved in 4 and 5 iterations respectively at Mach numbers 
0*9 with aft and fore effects, 

She CPU time varies from 4o sees to 1,2 minutes depending 
on number of iterations taken. 
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5 #3 SBwo-.ia.mensj.onaI IFlow 

Ihe ccmputations were carried out for NAC/i,- 00 l 2 airfoil 
and 6 percent thick parabolic arc airfoil. 

Pressure coefficient over NAGA~001X airfoil was 

obtained at Mach number 0 .VS in free air . She integration 

was carried out in the range « .28 < x < 1 ,28 and «9 .12 < v 

<9.12 and 39 x 6 rectangular elements were chosen with 

streamwise grid size A x ^ .04 convergence v/as achieved 

in 6 iterations talcing 2 minutes of GPU tine, She results 

31 

were compared with Lock’s finite difference numerical 
results, Ehey match v/ell near the peak, however, present 
is slightly lov/er negative at other points (see fig, 13), 

Results of 6 percent thick airfoil is given in fig .14, 

15 and 16, She integration v/as carried out in the range 0 ;< x 
< 1 "and w4,56 < y ^ 4,56, 25 x 6 rectangular elements 

were chosen with streamv/ise grid size A x = .04. Pig 14 
shows a comparison between iiie result obtained by Steger^^ 
and present method at Mach no, ,825 * Result for present 
method was found slightly lower negative than that of Steger’s. 
Oonvergenoe v/as achieved in 5 iterations taking 4o secs of 
GPU time. Pig. 15 shows the effect of solid wall on the 
pressure coefficient on the body at Mach number ,826, She 
wall was taken at (h/0) « 2,0 and ra,nge of integration was 
0 < X < 1 and ~2 <, y < 2 with 25 x 6 grid points. She 
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solid. wOill iirb 6 r enc6 givss iiigliGr nogabive pi?cssuJCG 
ooeff ioieni/ in "the mid region of the hody with 12 percent 
incscease at the mid point, ihe convergence v/as achieved in 
9 iterations with OPU time of 72 secs. ]?ig , 16 shows 
pressure distrihation in free oar for Mach numhers .7 and .8. 
[Dhe converged in 3 and 4 iterations respectively taking 
23 secs and 31 secs of GPU time . 

5 A Oonolusion 

An integral equation approach has been described for 
numerloally caJ-oulating inviscid subcritical transonic flov/ 
about slender bodies . It has also been demonstrated that 
the present technique can handle shock-free supercritical 
case. Ihe accuracy of the results obtained by the present 

technique Is almost of the same order as those obtained by 
relaxation procedures » 

Calculations with wall boundary conditions have shown 
applicability of integral equation technique to the study 
of wind tunnel wall interference , 

Further this technique can be applied to lifting 
proKLem, 
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Simplif ication of expression, 

E35)ression (3.11) is an integral over wall given Dy 


^ * Su 

^S.„ = “ C'f' I" dv d? 

-00 0 ^ ^ ^ 0= r 


W 


r p 00 

w 


w r w 

^ ^ WX) 0 P 


3,1/ p =5 


'w 


(A.l) 


No-w consider each integral over v ..separately with 4>=; -A- 

^ R 


S = [1(5 •0^+ ^ 2r r oos(e- v)J ^ 


w ■■ w 

Without loss of generality we tate 0=0. Henoe we obtain 

^ d V 

o R 




g(;t/g. k) 

% 


^(xIcTTTe,, + J')'^ 


where k = 


45 r 


w 


\ (x-O^Cw 

see 291.00 of refer'^nce jZ,^®I] ) 


(a . 2) 


Again., 


/ 4?p<iVi 


1 

271 

/ 

0 


oosv 

< dv 




1 

SIT 

C(^w 

% 

- H) / 

0 

av 

R 

% 

+ r / 

0 
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s;fC(Sw-f) 


1'‘(x^'7T(r^+5)2 [:(5-0®+(i„-Jf)^ 


~ k) — BC%^2 f k) — j 

'w . (?„ ^ ir 


(A.3) 


(see 291*01 and 291,06 of reference QES^]) 


Ifence expression (A.l) leads to j 


“ r.„ - r E( it/2, k) 


/ CK 


w 


w 




, P(7t/2, )£) - B(5t/S,]!:) - , ^ , 

+ 55 .,! ' " . . „ . ) fe >V 




^ P fw 
p % 


w 

E(7t/2,k) 


■] 


d5 




(x-C) + (r^ + r) 

A~II Differentiation of equation 19). 

3?lB differentiation of equation (3,19) w.r,t. x is 
straight foi^ard except the v/all integral term 


(A. 4) 


lo = / C(A+B) ?„)- 0^^; (5.y3« 

differentiating (A *5) w,r.t, x , we have 


On 


(A.5) 


& I 


/“c? ^ ^ 


a X 




-/ C? It f ST 


If 


Since 


9f (x-5) 


9x 

_-af(x-0 
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nere 


^ r 

df ^ ^ L 




® ^(x- 5)^+(5„+ r)' 


_ ^ k) 

r(x-o^ 

(see 7l0 *00 of reference [^28~] ) 


&nce 


^~+'“ 00 

[;0(A+B) / - / (A+B) dC 

^oo ^ 

[I(^C) U '^(x-5)^+(r +?)' 


« / (A+BfD) u(5, r ) d? 

i-iTi 


(A .6) 
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APPENDIX B 


So deduce expression 


Y/e ha-ve wall integral 

^ - I’ h 

^2 at wall 

= 6 H / (0 + / (^ tn ^r) 3 

(B*l) 

Since a-t c = + 

We haves 

00 i^ <50 

In = eC-/ + f ihr-^0r'> - 

^2 ^ ^c==-y^ 

(B.2) 

ipplyii^ wall boundary condition (3.26) v/e get s 

oo oo ^ ^ 

Iq « 6 C- 

■2 — 

= •-$/ {^0^(5» yyyr) t- 'I' ( 5 ~ ( C ""^yv^ Z] 

+ ?( E. y„) D? ( E =yw) - ( E 3 ^ (®-®) 

(since 0 and are symmetric about c -axis) 
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Now integrating by parts the second eapression in above 
integral and applying boundary condition (3,25) we get s 

oo 

= -s / { 5 iy„) { C *( ; = y^) ( c =-yw) 3 

- / p? c « = yw> « =-y„) 310 15 

(B.4) 
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iPPENDIX 0 


lo siiDplify ^iwSP expression (4,4) s 


x^+6. 


iwap 


/ (A+B+P) c. ? 


( 0 . 1 ) 


Integrating each term on R.H.S, separately: 




/ A S = 


r (r — r ) 


x^+6-- 

^ d5 




e ( V 2 , i %^) J ^;:^2 

1 


wi"bh a = (Xg— -f(rp-f?^) j b 4 Tp 


( 0 . 2 ) 


and k: 




= (h/a) 

^i 5=x, ^ (xg- xj_)^(rp +r^) 


Taking k^ constant over the small Interval Xj.-6^ l5<3CjL+6. 

so that E(it/2, k= ) Becomes constant over that interval, 

*i 

hence (0,2) gives 




/ A dE 




E(V2» I5.) rtan”^ itiiL. 

T % ^ (r « r. 


X.+6...- Xq 

Ic- » 

(j?w - ^p) 


1 X.- 6 .- x„ 

— 1 11b 




(0.3) 
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Similarly / B d? = _I 




2tc 




TTr 


In 


x,.~ 6,-- iXe + ^(x^-^^-Xg) (V Hp)^ 


1 ^1 S 


w^ '•P' 
(0.4) 


where l'(n;/2, k~ ) and Ii(n/2f k- ) are assumed constant over 
the small interval and 


- A 

^aq- ^ r v! 
^i-^i ^p 


B(n/2 


, )ta 




. ft ^ 


^+6; 


l^i 

(0.6) 


Oomhining (0.3), (0.4), (0.5) in (0.1), we ohtaan 


^iwap 


2tt " T. (?„- 


(V ?p) 


Ik^ ft 

^i"^i 


.-1 ^ 


l-k-t.,o l+ks ,e 


■tan~ - ^z:""" ^ +6 “ ^ ^ — riZi 

JU .1- 


E(V2. is,.) -E(n/2, f4+{?^+rp)^ _ Xg 


In 


Stc 


Y3 


x^+(v®y- % 


where » Xg - x^ + 6^ 

Xg == Xg - x^- 


(0.6) 
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ASSMDJX D 


®'i3SP 


2it 

^ f f f 

r 0 ; 


1 r c r , pdCdvflo 

*' '*'fx 


r ,-. q . 
3 ^3 


Integratioig R.H.S . w .r .t , 5 we get: 


r.+h_. X.+6. 


3 3 


^i3SP S- 


2 % 

1 j" f L P <iv dp 


^3 


0 


1 r 

BtT *' 

• ^r^3 


2 tz 


^ Z -A 


r X. 


/ 

O 


% 


L_(Xi + rj + p'^- 21^ 


2 - 22 - 
r X + r 4- p - 2r cos vj 
2 P 


3/2 


wl^ X ^ 4- 5 ^ 

Xe = Xg - % - 


Furttier 'we have % 


^ Ztt 


r + h^ 

3 

I * -/ J 
<13 ° 


(X^4- Tp + P^ “ 2fp P COS V) 


’ST^' 


(D.l) 


;s72 

oosv ) ' 

Pd V d P 

(D*2) 

d V d P 
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r .+h 

0 j ir dv p d P 

= J* 2X J — K ^ — 3/2 

„ o L A + rp + P^-firp Pcos vj 




r .+h. 

J J 4X ECtt/S, k) 

= -/ 


2 - .2 


dp 


X~+(rp+P) C2?+(rp -P)^D 


(see 291.04 of reference [l28 J 


where k 



ECTr/2, k) is complete elliptic integral of second kind. 
Further within ti^ accuracy of numerical computation we 
take 

P E(Tr/2 , k) 

'J 2 

X +(rp 4-f) 


constant over interval ^ 

as the point of constancy. We obtain 

4X r • ^ 

I =; ECir/S, k(X)) X 

2 


(X +(rp+rp 




vd- th P=: r 

j 


dp 

C (fp 



'^+(rj+ fp) 


r.+h - f 

E(Tr/2, k(X) t tan 

X 


- tan"^ LtlliiLfE “1 

X 


(D.3) 


P . 


( D«4) 
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where k(X) = ./ ^ ^.1 

Finally using CD.4) in (D.2) we obtain 

®ijsp “ k 

- 0(Xg) E(ft/g, k(Xg)) I(Xg):] 
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APPENDIX E 

I. Evaluation, of Ujj(x» r) for general parabolic body 
of revolution s 

The equation for profile of general parabftlic 
body of revolution is given by following two equations. 

The body for which maxiiauia thickness lies forward 
of its mid-point x 

R = 0 Cl (^*1) 

The body for which maximum thickness lies aft of its mid 
point : 

R=OE (e-E“) (E.8) 

max 

The following table gives the data for various 0 and n 
combinations with position of maximum thicknessj 


E “ 

(for max thick- 
ness) 

n 

R 

(the fornwla used) 

0.3 

1.71 

6.03 

(E.1) 

0.4 

2.S6 

3.39 

(E.l) 

0 .5 

4.00 

2 .00 

(B.l) or (E.2) 

0 s 

2.56 

5.59 

(E .2) 

0.7 

1.71 

6 .03 

(E.2) 


(see reference 11303 ) 
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The value (x, r) can he calculated using equations 

(4.2) and (4,3) with i being location of sting ( i es 1 

for no sting). However jbategral in equation (4,2) can be 
evaluated analytically for bodies having maximum thickness 
at mid-^)Oint. 

The profile is given by : 


4 R ( C - C ) 
max ' 


hence 

S( ? ) - 16 % ( ? - 5^)^ 

SUC) « 32it R|ax (25^ - S +5 ) 
S’*(0 « 3275 (6 5^- 6 5+1) 
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and iii) for x > jl 

I = 192,t eLoc ^ CS + I -13 


II . Fvaluation of u-j^(x, y) for parabolio airfoils ■ and 


HAQA OOXX airfoils ? 


For parabolic airfoil the upper airfoil profile 
is defined by 

(5) = 2t(5 - 5®) 

where is thickness ratio of airfoil. 

Putting it in equation (3,38) we get 

1 


y ) 


J_ / 2t (1- 2? ) d5 

It P 0 ^ 

2J _ r ( 25 - 1 ) In C(x-c )®+ 

Up*- “0 


2 «2 


1/2 


- 2 / Ui i;(S-5 ) + y 3 d; 


2T 


-2^V2 


^ + ' i;3^'+y"3 


1/2 


C (5 -x) In I (x-5)^ + y^] « 2(e -x) 


+ 2y tan”^ n ^ 

y 0 


2T 

It p 


-2. ^ 
X + 


( S' -x) ]ii 


+ 2 


2 ^ r tan”^ + tan"^ I 2 

y y 


( E.s 
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Similarly for NAOA OOXX airfoil v/e have s 

) ** Y? + 5 + Bg C + B0 C 

where Bjlj B2» Bg ana Bq are constants 


henoG 



(E.4) 

First we shall determine the y) surface of 

the body * Since airfoil is quite thin we can talje the 
linear Telocity perturbation on the surface as linear 
velocity perturbation caa the axis y #s o • So equation 


(3.37) becomes: 


0 ^) = If { s 


dY, (5) 


Using (E.2) we get 

oU = „t { i 


C— + B2+8B^5 + 

+ 4Bg 5®:]d (B.5) 


Integrating (E.3) term by term : 

J d5 7= 3^ f — R— where t^= K 

o 8(x-5) fC ^ o (x-. t'^) 

B. li- f S 

- 2zi — 

2 f- U fx 





U X 



X 


.1 SB A 
S __A_ d 5 
0 (S.5 ) 


884 / (-1 + — — ) d? 

« 2B4 1; -1 + S In 2 

1- X 


1 3B^C 


f :: 

o (x •*• 5) 


aK 


“SBg fZ - a *<)* -i-3 45 


1 


5c 

-2 _ 3? 


3B. 11- 0 - x + r IJi 


l-x 


and 


/ ^^8^^ 




d5« 4Bg / t- (5 + ^ )+ =?-r3^^ 

^ . (x -C) o X- C 


AH C"" ^ ^ j, 1vi ^ "1 

s» 4Brt t-** ST ** S' ’*' ^ + X -Ln .' '‘ "^ *''* J 


'8 


l-x 


Henoe (E,3) gives 4 


B, 


Uj^(x, 0 ^) ^ ^ Z ~ ^ 


1+fx 


2p 1-fx 


2B^+ 


(Bg + 2B. X + 3Bg 3c^+ 4Bg 2®) 
1- X ** o o 


^ + 5 ) - 4Bg ( ^ + 1+ r)3 

(E*6) 


1.x. 9 - 2 , 


Tbs abow solution for linear perturbation velocity on the 

plane of airfoil has sq.uare root and logrithinio singularities 

at the leading edge (x « 0) which gives inaccurate results 

near the' leading edge. It is modified by using following 
M 

formula * 

1 + ui,(5, 0+) 




FIG. 1 SLENDER BODY IN AXISSYMMETRIC FLOW. 



FIG. 2 THIN AIR FOIL IN TWO-DIMENSIONAL 
FLOW . 




WALL 


-ij 


3 LOCATION OF GRID POINT IN RECTANGULAR 
ELEMENTS ( axisymmetric flow). 



WALL 


JJMJljLlUJLL 
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X X X X X X I X 1 X I X 


X X X 

X— L--X- 


(xi.Vj) 

'thi 


4 LOCATION OF GRID POINT IN RECTANGULAR 
ELEMENTS (two-dimensional flow ) . 



i ^ BODY SURFACE ' 

, 1 I -.-i J J 1- ^ L 

0.3 0.2 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Axial distance , x 


6 DISTRIBUTION OF Cp FOR PARABOLIC ARC OF REVOL 
UTION (M„= 0.90, FR = 10 , free air case). 

o Experimental Present 

(Ref. 29) 











Pressure coefficient, 


“T -T 



Axial distance , x 


FIG. 7 DISTRIBUTION OF Cp FOR PARABOLIC ARC 
OF REVOLUTION (Moo= 0.9, FR = 10, wall inter 
fcrcnce case with porosity parameter =0.77 
and wall radius = 1.17) . 

Experimental (Ref. 29) 


O 


Present 






ON PARABOLiC ARC OF REVOLUTION (Mco= 0.9, FR = 10 




Free air case 



6 

0.77 


10 DISTRIBUTION OF Cp ON PARABOLIC ARC 

OF REVOLUTION (Mco=0.90, FR = 10, wall at 1.17) 


^ ^ ^ « 



1 DISTRIBUTION OF Cp ON PARABOLIC ARC OF REVOL- 
UTION IN FREE AIR AT SUPER-CRITICAL SHOCK 
FREE MACH NUMBERS (FR=10). 




12 DISTRIBUTION OF CpON GENERAL PARABOLIC BODY 
OF REVOLUTION IN FREE AIR (Ma,= 0.9, FR = 10) . 
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FIG. 14 DISTRIBUTION OF Cp ON 6% THICK PARABOLIC AIR 
FOIL AT ZERO INCIDENCE (M<„= 0.825) . 




FOIL AT ZERO INCIDENCE AT N 
IN FREE AIR BY PRESENT MET 
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'Jbo 


apfewdjlx k 


Ttilo (■‘HUGIMm t,;AliCuLiA'l’K{» EKPJ.'jfciUKE U'iSTKILUlX'IUN UVEH A UENERAXi 
EAKAlMllilC HlJDK UE KEVUljUTXDN iN EKt^ii; AIK UR IW FOKOUS wAUU 
W.J.t>Jl.l~Tll^'i^J);;l^ 


ii'JFil'r VARIAULES AKES" 

(M l i a i\l ( I tt * U !•’ F 1,1 1 N T » T A K B N U N B U I j U I, E IM G 1’ H B U 0 If 
•'iK^MUa.UF J.NTfe;RVAijS TAKEN IN KADIAB 'O-IRKCTION 
IjFkU'CJkB and AKT extent TAKEN J.N PERCENT UE BOUT IjENGTH 
MAXsMAX IM'IM NUS, UF ITERATlUNfj ALiBOwEP 
J.tta:U,EUK I'lU STXNU 
1,'VlTH BTING 
XaTsijUCATlUN UF STING 
KHSKAUIUS OF KTING 
Ilj=i.),E()i< free AAR CASE 

1 ,i.'Ukwaijiij antbhference case 

FKs fineness ratio OF THE BUOY 
NsFREE stream HACH ND 
AEUaEUlJATAUN NO. USED (APPENDIX BJ 
CU,PU ARE CUNSTANTS FUR THE BUOY USED 

rnkkaijius of wind tunned WADL 

wpapURUvSlTY PARAMETER OF WADD 


THIS iS MAIN PRUGKAM 
READ M 
INTEGER S,P 

CUMMI)N/lNT/S,F,X(t>U ,R(10J ,H(10) ,NX,NK,FC50f 10) ,RB(bU) , AB(5U,10) 
CUMMON/vSE.l/ l)E,|JYl ,UYi!, IST.NP 
CUMMUn/SBTI/DEIiX , 1U.BL4,IS 

CUMM0N/SKT2/RMAXf Y/R5 
COMMCIN/EF/EfFE 

CiJMMUN/XW/Aw(40,4U,10) .UU4U, 10) ,WP,RW 
COMmUM/BUS/CO^FU, lEU 

DIMENSION ODLCSUflO) .UDC 50 , AO) , UND ( bO , lU ) , 0 ( 50 , 1 0 } 

1 ,A 1 N( 50 , JO j;,EU 0 U,FE(l 0 i),DH( 40 ) ,CPDC 40 ),CPNUC 40 ) 

T ^ Pii y 5 0 

FORMATC lOAir 'TYPE «B NR liP MAX IS AST RS IL FR H IEU CO PQ'’) 

ACCEPT ;NH/DP,MAXf IS,XST,RS ,ID,FR, lEQ^CUrFO 



"* 'I 


iKtiri.MK.o) TypK ysi 
F UK f 1 A 1 ,' c 1 0 X , " 'C < F !<; K W W F *’ ) 

I {<■ t I, ti , i'iF , U ) ACCFt'T , K W , WF 
P K ,1. X V J. U , N W , (') K , jj P , M A X 
PKINT HXUpl;:),xaT,K5 
PtUin' (KU), IbjFWfi'l 
9 K .1 X T « .U , X fJQ / CU , P CJ 
X F I I b , X F . (J ) P K I NT 8 4 0 , R iv , W P 
FUKilATCaX, 'KWs'^Fb.J, JX, "WPa' ,F13, J) 

FUH l'1 '' I' ( / £> X , »■ NM= " , 1 2 , bX , ' NK= *• , X 2 , 6 X , ' LP= ' , 1 2 , 3X , ' MAXa ' , 1 2 ) 
FUKMAT(bX, ' iSaSH ,bX, '’xaT=:%Ff),3,3X, 'Raa'#FV»l)J 
FURMAT(bX,'i{,s',Xl ,7X,"FR=',Fb,2,4X, 'MS' /Fb.i] 

F U K M A 1' ( b X , ' I E U s » , 1 1 , 7 X , ' CU = ' , F (1 , 3 » 2 X , ' P □ a ' , F 6 , 3 ) 

IFl I.S,NE,U} PKINX yjO 
IFC XL.Fil'U.U) PRINT 920 
1, F ( .1, 1,1 , , U ) P K X N T 9 3 0 

FUKMAr(/lUX, 'rtUuy WXTH HTXNU'j 
FUKMA'I’C/lUX, 'FR fE AiR CASbi'J 
FURMATC/lUX, 'WAbb XNTfcRFEKENCK CASE') 

EKKUK=O,0Ul 

YsXttT 

RMAXsU,b/FR 
(3al,4 

91 FuKMATClOX,F7,3,l2Fy,4) 

UEFjXsI ./FbOATCNMj 
Wi3 = l .-MI'M 

CUWasKH'l, ) ♦MtM/MB 
iVF = NMVL,P/lUU 

X U .) = U . S-FbOAT t NP ) .) >«REbX 

iHaNB" 

XFi; LS,NK»U) I»sx^KIjUAX(WB)+U,bl 
IVXSi,B + 2=I‘WP 
XFCJS.CJT.UJ NXaNX + 2 
jaX'slM + NP 
XUSTjaO.U 
LMj 10 la2f<'JX 

1 F ( X a , EW .1) XX iaT+ 1 J ax ( J. ST ) 
xxiisxx i"‘i)+ofe;bx 

10 cuntiniIe 

XINX+1)=XXNX) + DF)jX 
XFCXa.EW,U) GU TU 11 
aTaXtlST)+U,£>VDFGX 
UEsY-fiT 

0XlaXUElUX+0F)/2* 

UY2al)EbX-UK 

XdStjaXC.lST-D+U.b^CPEbX + DYir 
X(X5T + 1 JaX(XSi’) + f)Yl 
X(XaTt2jaX(x5'X'+U+U,!>=|'ClJ3fl + UY2J 


ybl 


b4U 
H I U 
tt2 0 
«iU 
« i 1 


y 1 U 
y2u 
y.ju 



i 


(NXaiMX»l 

ji 1 (juM I.’ I KHJh; 

Hin=:Kl’<W\X 

K(n=H(l) 

LUj XO vJ = 2,WK 
MI..J 

l< ( t J ) “ '< C ^ J “ 1 ) + H I i.I "• 1 J +■ H ( J ) 

iju Ctji’jx.iwut'; 

t K l J. 1 . 1 , r»lb; t 0 J H U^K) = ( KW"( tu MR "1 J + H ( NR-l J ) ) ’t'U * h 

IK C .1.1) , l^lb;, 0 ) R^NK}=HINh*» 1 J+HlNH-l ) +H INK) 
un .51 .1 = 1. , NX ■ • ■ 

X.L = XLl) 

C A hU W U KF U 1 f K vj , DIU , t) , 0 ) 

KflCi)=KJ 
l)H( 1 )=Uin 
uij1=U):kxi;) 
ujiLU , nsui-ii 
ujtC. i , j ) =111)1 vcuNa 
UK K U=’Ul,.iC,i,l ) 

KUKl-iiKliD^UiaK) 

i.>U ,5 0 UK';i,WK 
KvJ = KCJ) 

CaLiIj l)XNEAKUl,HJ,ljLl,CKj) 

UIjIjU# JlaUlii 
IJ.U( i! ,vl l^lJUKCUN.a 
UK KU)=U1)(1,U) 

KCl,U,)K:U.l(i,J)'KUKI,vU 

,5 0 KJWTlNIJh; 

3l CUN'IMNIJE 

AIWIjSU.U 
A INTO a 0,0 
ini bU 1’1'=1,MAX 
lJU ^11 8 = 1, NX 
OU YA Fs1,NK 
L‘5\bU IN'VKUI A.lNTU) 

iKC (XU,NE«0J ,ANU» IK’.KQ.I ) ) CAljL WAIjIj 
XkC CaLIj XNTWACAIWIj) 

A.'lU2«"’Ar5(a,F)/3 t l<lisyk!l>b4 

UlLS,FjsiULCaf F)"'AiNtG-ANll2=5'K(K,F)+AJLWlj)/(l*"ANU2’l'UKS,P) ) 

AJ.N la, F i=A.tNi’ci)"A.iwit 

YA CUNfXNUK 

CUN'KNtJK 
UU bl 1=1, NX 
UU ‘jO J = i,NK 

XKCA»a(U(,l,U)«l:}KX, J3 ) ,«T,KHKUR) (JO TU lOU 
1>0 CUN'KNUI!; 

bl cOntxnuk - 

FHINT y2, IT 


I 



4 


7, I'' * ) l< i-l A i; t HI X , .1 ^ ' T H I T E H A T X U N ' ) 

UU I.= l,iMx 
MU :>.i ^) = l, 

b.i L'UWTXMUK 

UKUaU.O 
UhXslIiKX) 

IJKH:s|)KX'i=DKl 

CMti 1 1 ) =«2 , C 1. , X ) -DKK 
Cl-'iN l.i 1. I, ) =->{ , VUIM l,'i 1 1 , X J -DRK 
UUNI'l I’JUI:: 

I'Kl NT 9i, iKXJi , 

y i FUR M A i ; (/ 2 7K ,' K = HUP'/ «i) HF ' , BX , ' K = ' , Fb , i , 1 UX , " Ra ' , Pb , 3 , H>X , 

1 ' K = ' , Fo * 3 , 1 OX , , Fb . 3 , lOX , •’Ra ' , Fb , J / 1 4X , ' X ' , 7X , 'UIj " , bX , 

, /X , '■ur.i' ,bX, ''UNI j' , 7X, HJlj' , foX, 'UNL' J 
FMI.I'IX' 9.1, UU. J , lUMj,((:X,vn ,lliMl.iC.I., JJ ,J = i,.NKX ,1 = 1,NXJ 
IMT 94 

94 Fi)Hi'lAI,'{/luX, 'tJCIKFF, OF FKKSb, l.liSl', UlM tUJUV SURFACE'/ 1 4X , •'X ' , 

HJX , 'UF-Xj J.WFAK' , 6X, "Ct’-NUi'l ij.lNfcWXR' , 4X , 'RFt) , Fb'Ki', VFL, U ^ J 
FK.tNT 9b, UXti),CFi;itl^,LPNHIJ,IJ(I,Xj),Xal,NX) 
y b M AT ( X 0 X , F7 . 3 , bX , X' 9 . 4 , 9 X , F9 , 4 , 7 X , F9 , 4 ) 

iSTOF ' 

XOO UU b2 X=X,NA 

UU bX ,«i< 

Ml. (.1., J)=UU 

F U , vJ J ~U U 1 , M X >I‘U .t t .1. , J ) 
t>l XaiNTJiMUK 

t>/. L'UIVTXNUF 

OU CUIMT.l.NUF 

STIIF 
FNI) 

x; 

C 

U THXkS iiUBKfJljtxiNE CAJjUUPATES 'XHF HOi>Y' RADXUS ANU ITS UEKXVATXVES 
'(: at DESIRED IjUCA'nUtViS 

SUBRlJUTANfe: SURFU,H,DK,DI.)K) 

UUMMlUV/SET’-i/RMAX , X ,KS 

■ CUMMUN/BUS/CfF, lEU 

ZiKRaO.U 

kaZl'iK 

DRk/^FR 

liURsiiEH 

IF [X.liT.ZKRj RETURN 


RaRH 


j,F(X,UT,Y} 

KETUKN 

Ch!=c»kmax 


tt 

X 




cv n n 




ill 


l: 

(j 

L 

c 

t: 


c 

t: 

t; 

C 

i; 


i: 

i: 

C 


ll-'C 1. F;;u.|';q, 11 ti=!:U"X 
Al.J=:u,U 

1 . (* ( a f f I 9 0 a F) 0 1 , ) c^ll IH) 111 
= U’”-! a J 

1 a-AlSJ 

UI>H=«L'R^=pvip-.i , HAB/B 
UKSCW)? ( i .-F^AlU’^t-l « )>P* 1 ISU 
ki';ti,iW(\( 

l!,Wl) 


liUi'i t Ui'IC riUN UAl,iClUjATfc;,S OtKiVAlMVl!. UF -BUPif tJKOi>S"i>I:;CTlUNA!i 
AKKA AT AXlAh LOCATE UN X 
IfUWCTilhM USlXl 
UAbU BUKFtX, K,1.)K,UDK) 

UMkk/;; , vHK 
KFTUHN 
h; i'i 0 


THltt (■'liwCTXUN CAriUUliA'i’Ki) {jECJOND UKRlVA'IlVE UK BODY CKUSS" 

i>«;(:.i'j:i;wAi, akea at X 

ruNCTiuiM iJiiatA) 

CaIiIi aUKFlX,B,),)K,UllK) 

D t ) a = U a b ( U K V UK + K > 1 ' U U K ) 

KH.TUKN 

END 


i'KXW Ifj AUXifiiAKl* KIINCTIUN IN CALUU i-A TlNG (JH ( X ) ' 

FuI'JC’T.IUN FC-TCX,fc;) 

FCTaO.U 

UEKsAlWUX-El 

IFIDU', IjT.UfUUOl) UEXURN 

i'--Ut)«U J-UUrtCh:) , 

FCTai'VUIK 

KETUKN 

b,w|.) 


if-lTW fUwC'ilUN CAUCUUATEb biNEflK PERTUKBATlUN VELOCITY OVER 
Bum auKFACh. ax axial lucaxion x 
^ UNGTltlN UB EX J 
CuMi’JlUN/yKTl/UEEiX, .I.B,BB,IS 



tin ^ 


U t.J M .VI i I f'( / » (-; '^ / H M A X , Y , K H 

fcil/M = (),() 

It. J = 0 , 0 

uu li) .rsi,jH 

l‘ .i< = t;;n U .SY'Ubilj 

h i M p s (•• c T ( X , h; 1 j fif c i c X , ii; j ♦ 4 . 0 + 1 ' ' c j; ( X , w 3 j 

fc<U|V|3,'jli)'/l+'SlMh:^ 

lu L'lJWVI Mill'; 

fc'UM=i.)i!;):i/6,*rtUiv) 

Jlt' ( .l.b.tCU, 0 ) Gu TO luo 

tttiM=!SijM+ ( Y J / ^ , CKCT ( X , Y } +KCT I X , Ei J J 

iUU GtlW'l/lMUK 

hliEKKUSiyj/AHHCX-Y ) 

LALlt aUlU''U,H,|>K,I.Il)K) 

JfCX.GT.O.U , AW|>,X, fjT, Y J Sl..l':f-BliKl'' + DUS U ) ♦ AliOG t UlU^K’I'K/ C 4 , C Y-X ) ) ) 

KIli'l'UKW 

h,Mf.) 

C 

C 

C 11116 aUHIlUllTlME CAliCUl.ATEy l.lNHlAR RERii.lKBATiUN VEJjOClTY 

C AWAY f'RiJi'i I'UlOY «UKFA(>; 

ttUHROUTltMti; JjXNEAH«.X,R,Uli,l.’l’lj) 

CUM ivi I.) M / s I*; T 1 / II h; L X , j. n , n h , j. s 
CUMMUW/»El/ l)E,UYJ,DyX, 

1 1 1!) i.t 1 1 Kt li X / 9 

HKsBMI'RVR 

WiiMKlI.U 
i ll I 4Ul) J!3l,l2 
l;;aE + Ul’;i.) 

I'Rsil , /SURtU-)*M + HKl 
U t.i s U / S U R T ( 1 U + 1) E ).j J ’I' C I ) + 1 ) It: I J ) •«• H K 1 
SUMs6LIM+1)UC EH (BR«U«) 

4 011 CUimMlJK 

XEt US,NH;,1) GU XU 100 

a I j M s: a LI M + n H t E ) 'I' ( 1 , / ■ j '■>■' t ^ ^ ® ^ ^ fi ^ ^ ^ ^ ^ 

100 CliNtiNUE 

0CS6UM 
ijRLis-ii'.’l'Ul.i 



n n n n 


kktunn 

KNI) 


■A'Hia aUHHUUTlNE EVAliUATES THE ElKUD INTEGRAL ' ' 
tj|JHR[jUTj.NE INTEU'CAlMTli) 

INTEGER S,P 
REAL KK 

CUMMUN/i|«JT/5,H,XC!>l) ,K(10) ,HC 10) ,NX,NK,E(!5g, 10) , RH ( bO ) , AB C 50 , 1 0 ) 
CUMl1tli\|/8El/ DE,UYl,!)y2y ■ XaT,NR ' 

CUMfiaw/8gTJ,/l)Er.iX,,lH,Bli,IS 
L’l.iMwnN/EF/E,FE 

iM«EN8ltlN A(.bl,lO) ,HRtl0),BH(10) ,E(lOl),FECiyi) 

NAlsNX+l 
tisSvSURTCHR) 

BKRssHYRtp j 

J:F(F.EQ,1J HHRaH'KKlHS) 

DU y vJai,NR ' 

«H(U)aH’»'iKvJ) 

HK(vJ)=B’I‘rCv)) 

y CUMTINIIE ■ 

SUMaO^O 
UU 10 1=1, NXl 

1)eLTkuEI.iX 

IFClS.NE,!) GU TU H 

iFC:i,EU,I8T,UK,I,EU.CiST+l)) UELTsDYl 
Xi?’(,UEU. (iST+^i)) DELTSDY2 
a CUMTiNUE . ‘ 

X«IaXtS)«X(i) + 0,5^I)ELT 
iF(l,E«,i) GO* TU 20 
BKlsB^RBC 1"1) 

KRLaBKl+»KF , 

KKPaKFL^RFL 
XK aSyRT(XSl!1‘XSl + KRH) 

XKMa«URT(X8H»XSM+KRF) 

AK=4,U»BRF*HRJ, 

KK»AK/ CXR*XK) 

CALL ,LNTPUL(KK,E,EK) 

mxsak/cxrm^xkm) 

Call iNTPULlKK,E,EKM) 
kMI'J=BKl-0KF ‘ 

BHKa2,0»BHCl)"BKl 

kRNsRMW»C«Hk+KMNJ 

AiSa8KX^(EK/Xk'»ATANCBHR’('X5l/USI*XSl + RRN))-EKH/XRM»ATAN(BHR*XSM 

l/CXS«*X8M+KKNj)) 

AH(X«i, i)aAiB 

SUMaSUM+P t 1 »a ,1) »AB ( I-l 7 1 J 



n n r: ' ci^ cv cv ci r: 


" H 


•^{) CU^IMWIJK 
X ti M 3 X H I 
UU AX 

KHI,tlS = BK(iJ) + HKP 
A K = tt Q K T I X » 1 ’i' X a I + R F I j U S 't' R P ILi U S J 
M^ = '1,0»MRP’^BKUJJ/UK»XR J 
CAr.ili .!.NTPUl,rlKR,l5;,KK) 

KMNaahRUU~WRF 

At :i, ,iJ JKHR UU /XR*l!,’K»>'CATAN( I RmNH+BIU J J ) /XSI J -A’l'AN ( ( KHnS"BH ( J ) ) /XSI J ) 
iKCl.fiU.lj «0 TU 21 ' - 

AH C i ".I , J ) aA U , J J "A ( 1"1 , U ) 
bUMaauM + Ftl-l , ,iU 

AX CUi'J-.ilWUl!; ' . 

lu CUNTlWUl?: 

Ail'JTGs.b/J, x4l5V2bb4H<a(JM 

KtiTURM 

Ki'll) 


'i'Hia MUHRIJUTXnK iNTKRFOtiATKS ijINEARUY A FUNCTION BETWEEN TWO 

PUIlVTH ■ 

BUBKOUTINE iNTPUtiUKjFCT, VAl)} 

DIMKNSXUN FCTUlUlj 
C K K s 1 0 (I » ’t* A K 
M=CKK+ 1,0 
OlFaCAK-FUOATCMJ 

VALaFCTtMJ + tFCT(M + l)-FCT(«) J^IJIF 

RETURN 

ENO 


» ♦))(;♦! .K ^ )j! )f ^ If 3<< !« JH Ht f f ♦ 

THIS aUHRlHJTiNE EVALUATES CUEFFiCiKNTS FOR WAlUli ' INTEGKAIj 
SU lJKflUTlNE WALlJ 

iNTi'JGER S,F • 

ki*;au RK»KX,M1 

CUMM»l\i/XWVAW(4U,40,101,Ul(4U,10JI , Wp,RW 
CUMMON/BBTl/OEtiX,XB/BB, IS 
CUMMON/aEl/ DE,UYlrOy 2 / IST#NP 
CuMMUN/KF/K# FE 

CUMWUN/iNT/a,P,X(!jl) ,K( 10 J , H ( 10 ) / NX# NR , F C bO , 1 0 ) #RB{ 50 J , A6(5U# 10 ) 
DiMENSlUN Ell0U#FE(10l) 

BbSUHTCUB) 

BRWBOfRW ■ 

BHPatJ’I'RtPr 

IFCP.EO.l JB kPbB»RB(S) 

bposbrw+bkp ■ ^ 



« y 


b 


lu 


12 


XaisX <. S J -X 1 1 ) +0 . 5 >i'OK/,X 
y.I. = «QHTtXSi=PXai+BFL¥BFL) 
Xi<'(hHP,UT, 0 .UU 0 l) GU TO IJJ 
BfciOsSUKTChKW’titvHF) 
M.>lNaBKW-BKF 

fc 5 K F . 1 . = 1 1 R W / , / 3 » 1 'U b y i» 6 5 ) 4 

TX = ftT>\MV-X{jji»M/BMN) 
rtliX=Al.UU(-Xtti + yi) 

C‘ljX = AbOtj( U 1 ,+M) ) 

uu 10 1 = 1 ; NX ■ ■ 
ii = ,l + l 


ije;LT=uKi.iX 

XeCia.MK.lJ tiu TO 8 

I r ( J. 1 , I!;m , J. «T , UR , i I , KQ , ( XvST+ i ) ) DKliT^Oy I 

ir c 1 1 , t:u , i JLBT+2 ) ) l)b:LT=l>y2 ' 

CllMTlNUt; 

x«.i:i=x CH)-xui )+o,5VDi;;iiT 
y X 1 Si B U R T ( X S i 1 » X i> J . 1 + B P i I ♦ B P | j ) 
K.u,= 2 ,»bsu/yii 


'.lM.iaATAN{-XSIi>J‘Ml/BHN) 

Ai;;a=AijU(u*xsii+yxi) 

UblissAIiUGC C l .-Kll y/(l« + Kll ) J 
K K = 4 ♦ R W * liK P / c 1 X L S ) - X C 1 ) ) » 2 •!• B P U » B P L J 

CALiIj iNTP(Jlj(,KK,l!;,t;i) 

CaIjLi XI'iTPOU(KX,i''f::f Fi) 

ifC tKX-.gyOUJ .GT.O.OUOOUIJ l'l = 0,5>l'AlaQ(»(l6,0/Cl ,-KX) J 

AWU. ,»,P)«Ul/B!30»C2.H‘(TU"T.l)-WP/0»(CLli»C|j3:) ) + 2iO/'BRW'f=(FI-E;i) 
1=»‘IA011*''A14 ) I^^HRPJ: 

XBlaXBll 


■yl=y.ii 

KjalVll 

TiBTJi 

COXaCLiIl 

AIjiaAiili 

CUNTiNUb; 

(.ill TO 21 
eUNTlNUbi: 

DU 20 t«l ,NX 
1 1=1+1 

DKI(T=UE1jX 

Ib'ClS.Ni!:. 11 GO TO 7 

JFU.1 .EU« IST.UR.Xl.EQ, ClSX+1) J UEDTaOn 
XF(ii«EU,UST+2)) 1>EM'»DY2 
CUNTINUK 


yllaauHTCXSil»XSXl+BRW^«BRW) ^ i / v t i *ysx/ VT 1 

■ AW c i , s.p) »o ,s+ iBRw/B=<'WP»c 1 « / 7 i )«xsii/y Il+X&I/XX) 



- 10 - 


^ fc* I = X M 1 1 
X 1 = M 1 1 

•/JU Cui'J'l.'itMUH; 

2 1 Cun Tit'll I hi 

KhlTlJKN 

lilNiJ 

c 

c ■ 

C; TMIH JiUBRUUTiNK CAIjCUIjATES tHhl WAIjIj ' I NTKGUAIj 

auUKOUTlNK iN'l’WAtAlWij) 

INXyCWH S,|J 

CUM.^UiVJ,iMT/a,P,X(bl) ,KC10) ,H(aOj^WX,NK,KCby,iO) ,Rli(bOj ,ABC5l),103 
CUMf'1UN/lW/AW(40, 40, lU J ,UH40,iU), wp,RW 

Bl.l.v| = U,0 ■ 

Dli 10 I a I, NX 

^UM=aiJf<1 + l.)l t J. ,NHJ*AW(I,S,P) 

HI CUNT l.wUtt 

A.LV'UiaailM 
KhlTUKN 
WIVI) 

c 

c 

C hUljljlIWlWG AKP: THK VAtjOeti OF CUMPLiSTt: tlLOlPTICAL iNTKGKAliS OF 

C ShlCUNU AiMO FIRST MNDS 

BOUCK data 

DiMFNSIUN PU1UU»FE(IUU 
ciiMMON/iir/t'^Fii; ■ 

DATA hj/i,b/U7y6,l«5b6H6:i,l, 561^913, 1.5b8948,l,?)549b9,l,5S0y73, 

11 ,a4fa9b:4; 1, 54i!y3b, 1.638893 ;i,b34a33, i,6i07SB»i,526b65,i,b22I>b5, 
2 1 , bl 841^8, i » 5i4i 84, l.b 1 011^2,1,605942, 1.601743,1,497 626, 1,493290, 
3 1,4890. 16,1,4847 61, 1,480466, 1,4/6152, 1, 471817, 1,467 462, 1,463086, 
41,458688,1.464269,1,449827,1,445363,1,440876,1,436366,1,431832, 
51.42/274,1.422691,1,418083,1,413450,1,408791,1,404105,1,399392, 
61.394652,1,389883,1,386086,1,380259,1,376402,1.370515,1,365596, 
7 1,360645, 1,365661, 1,360644, 1,345692, 1,340505, 1,335382, 1.330222, 
81,325024,1,319/88,1,314511 ,1,309192,1,363832,1.298428,1,2929/9, 
9:iL.28/484,l,28l9 42,X,276360,l,27O/O7,i,265Olf, 1,259263, 1,253458, 
1 1,247595,1,241671,1,236684, 1,229632,1,223512,1,217321 , 1,211056, 
21. 204714, 1,198290, 1,191781, 1,1851 83, 1,178490, 1.171697, 1,164/98, 
3 1,15/787, 1,1 60656, 1,143396, I, 135998, 1,128451, 1,120741 ,1,112066, 
41, 10477 6, 1,09647 8 , 1,087937, 1,079121,1,069986, 1,060474, 1.050502, 
bl.O39947,l.O28595,l,Ol6994,l.O0QOO0/, 

IFF/ 1,570 /96 ,i ,574746 , lib /SVAO , I » 6827 80 , 1 ,586868 , 1 ,591 003, 

2i ,695188(11 ,599423,1, 6037 10,1, 6O'0O49, 1,612441 ,1,616889,1,621393, 
31 , 626955 , 1,630575,1,635257 , 1. 640000,1 , 644806 ,1,649678 , 1 .654617 , 

41 1 6596'f4 , 1, 6647 0 1,1 , 669850,1 i 6 /5073 , 1, 680373 , 1,6857 50 ,1 , 691208 , 



hi .t;>yt>74y, ] 
fii , /3wah5, : 
71 

M 1 .H .11 /4yi »] 
y 1 , MyH9'^h , ] 

1 1 . y/i7H3 , ] 
Zi?,ybo8u;^,; 
3/5. i 7 4Hi57,; 
<»'/5,.J3l409,; 

fi3, IbhHVb,, 
(i;''!)'; 


1 ' 

. /443bl,l, /hu /b4,1.7b7i6«, I, 7b3byb,l, 77^047, X.v/Vbao, 
. /ytbbO,!, /g«y 18,1, y0632a, I, «U8a4, 1,821 b93,l ,«Zy460, 
.84Jby4,l,a54075,l,«b2B4i, 1 ,«7i4UU,l,Ha(7i6l ,1.889b3i, 
«9u8b47 ,i,yiu4io, 1,928526, i .yiayoa ,1,949560, 1,960521, 
.9833 7 1,1, 995303, 2. 007596, 2, 020279,2, OiSiby, 2, 046894, 
.0/53 63,2,0903/3,2, 105948,2,122132,2, 138970,2,156516, 
,1939 /O, 2, 2 14022, 2, 235068, 2, 257205,2,280549 ,2,305232, 
,359264,2,389016,2,420933,2,455338 , 2,492635 , 2 , 533335 , 
,62/ / 73, 2,683551 ,2.747073,2,820752,2,908337,3,016112, 
,354141,3,695367,999999,9/ 



X 


c; 

c 

c 

c 

C 

i; 

C 

C 

C 

C 

u 

C 

'C 

C 

i; 

C 

c 

C 

C 

C 

C 

C 


yi 

y&y 

y&i 

wio 

B'^U 

BiU 

yiu 


i'Hi,'j k'HUGKAM CALiCuljATlfJS PKtlbSiJKfc; i)lbTKittUTit)N' UVEk ft BAKAttUblC 
AkC AAHJ'OlJi iN IfKfcE AIR UK JlN PURUUS rtAbb WIND TUNNKb 


INPUT VARXAbliEaj- ‘ ‘ ■ ' 

WUSWUfj.UF PUInTB taken UN AJLRfaib CMOKU LENGTH 

NKSNUto.uP PUiNTtt taken iW LATERAL Oi’KECTiUlV 

l/PrEUkh; ANU apt EATEN'T TAKEN TN % OF AlKFUlL LENGTH 

MAX=;MAX.I.MUM NUS.UE ITERATJLUnS' AbLywEl) 

TsiTaiUKNE&B KftTJlU UE fti.KFUlL 
ivisT'KEE stream MACH NtJ, 

J.|iS!U,EUK EHblE AIK CASE 

1 ,H'UK WALL INTEKEEKBNCE CAKE 
)[W5!HAI.jF height UK THE WXNU TUNNEL 
WPapukUttIXY parameter QK WALL 

f ^ )f if >)( V 4 : V V >!> >K V V >(( V 4^ 'F V f V >K * « V r r f « V V « ¥ V * V « V )« y ’ll ¥ If 


'J.‘tU.;:> iS THE MAIN PR UGH am 
Real M 
iNTEGEK S,fc> 

C()iLnUN/Ul/S,P,NX,NR 
euMivHiN/u2/UlC^0,lU j ,Wl4y,40,lU) 

" CuMMUN/UJ/Xi4U) ,YClO) ,HtlU)/WP,»,UELX,5fW 
CUM RUN / U4/ K 1 4 O , 1 0 J , A ( 40 1 1 0 j , Y«PU { 40 J , K' BMN ( 40 ) 

Cumrun/ty/X' ’ ' 

CuRMUN/UU/ULLUU#lOJ,UNL(4y/iy),M,G 
DcMENKIUN U14U,1UJ , UBC40,1UJ ,AAUU) 

Kuhrai (, iOX,K7,3,l2Ky.4)' 

EKRUK~y»OUl 
T X t-* i?i y 5 u 

KUKRAtCiyX/'TYPE NB,NK,LP,HAX,X,rtriL^ ) 

ACCEPT 4‘,NB;NKrLP»MAX,T,M,IL 
XK(,.i L,NE»y) TiiPE ybi 
KiiKMAT(lUXf'^Ti(PE ' 

TKlXL.N^.yrAeCEPT *,YW;WP 
PkIwT 8iy»N0»NR,LPrMAX 
PHXNT 820,IL,T,H 

IKUL.Nf.U) PRINT 0iO,YW,WP ^ .««.v 

KURRAT(/t>X, •NbB*’ #l2,6Xi 'NKs'/lZrBK# ^P-/ / If » 3X r ^ MAX 

KUKRA'i (bXif *’YWa*'f Kbi3# 3X, 'ftp*'7f'b*3) 

XKCiL.EatU) ■PRINT 'yio'^ 

XKC1L,«E.U) print 920 ^ 

KUKRAT ( / I OX / ' T-KEE AIR CASE ) 





*n 


2 ’- 


'u/\i,h Ki-’ir:f(C’«: cArtf','") 

WK~.MS • - 

i'i K .I = i 
f4i<2 = 2 

0 = 1, 

L)l?;La= ! ,/K»jOAXlNHj 
IrlaWljK'H 1 , j 
uuwasui+i . ) (b’I'B) 

Wt'SNH’f'IilVlUO 

Xai=5U;.5-F(jUATlNFj )»L)E|jX 
rUKWB + X'fNH 
uu lu .1 = 2 » NX 
X u. i»Au-i 3 +uk;IjA 

10 ouNTiin’ih'; 

!«* t J.N, wf!;.oj Gu TO n 

u.1. l=lfMAX 
uiii = r.^YMAA 
‘HU1«K{1) 

Hl2j=i3.»YMAA 
'UU 20 J=J»NK 
Ht J)=H(J-1 )+5,»XMAX 
U J jail J-i ) + Ht J"1 )+HC»J) 

20 CuMTiiMUt; 

OU i'C! 22 

11 1>W = U, !>=f iw/ CNR-n 
HllJ=UH/2.* 

U1)»M(1) 

H(2)aDH 

Xl2ja2,^DH 

UU 21 ' J«3fNK 
ii'lJJPIlCU-1 J 

ilGjaX ij-U + HlJ-lJ+HCU) 

21 continue 
BcNR]=MUJ 

i (iVK)=YiWH-a 3+HlNK-l) + HtNK) 

22 CUNTiNUt • • • ’ ■ * 

UU il 1*1 rWA 

X1=AU) 

XfcjPli(i)=YX (All 
■ YJsO,O 00U01 
UbCl, 1 )«CUNb!fUtitXl,YJl 
Fit, 1 j=ufvU» i)*uijiivi j 
ulc.l, IJaUb tl; !) 

Du 30 0«2 I NK 

UTi i *’>")=CUNS*ULiXl ,YJ) 

f Jr 



m 3 ** 


U1 (I,vnauii( 

JO cutiTiwuf-: 

ji cow IT wilt: 

y ( 1 )»o.{) 

AiWIjSU . U 

Dtl 00 X'l’sl,MAX , 

DO (11, Ss:l,NX 
DU i\2 Psi,l\!Y 

CAlili VW'fKGCMMTC;) 

;i;p'(,r.t'.ti:D.i,AND,i.L,MK,oi cal,i. wai.l 
rp(i:r,.i\ji‘;,a) caIiL .iwtwacaiwi,) 

AACP)=AXNTCl + AIWb 
Al\!ir,^s:A(«,P) 

I.IC«,P)3(llBCS,P) + AINTG-ANU2*F(S,P)+AlW[j)/(l *-"ANU2*Ul(S,P)) 
42 CUWT.rwUK 

4t CniMT.XWUtl 

DU T.5fl;NX 
DU L)U Jnl,Ny 

lF(ABvS((j(I, ,GT,ERRPR) Rn TO 100 

bO CUMTIWIJK 

bt cuaitiwiip: 

PHIMT 92, XT 

92 T'’UKMAT(10X,I2, 'TH ITERATION ' ) 

PRINT 9i, CYCJ),J5l, NR) 

y;3 F'ORMAT{/2 3X,'ya'',F6.3,lOX, '’Ya'’,F6,3,iOX, ' Y**' ,F6 * 3 , lOX , 

1 *'Y = ',P6.;J,10X,TY='‘,F6.3,10X,'ys',F6,3/14X,'XS7X, 'Uli', 
26X, '’(INI.%7X,'UL',6X, '’UNb',7X,'Uh',6X, 'UNl.'',7X,'’Ul.'’, 

3 OX , ' IINL" , 7X , " UL' , 6X , ''UNL ' , 7X . 'UL % 6X , UNI.'' 3 
DO b? laJ.NX 
DU bJ 0=1 /NY 
Di\i|.CT ,inaUC j:/0)/CQNS 
lO'I.CT / J)siUMI / vD/CDNvS 
5 3 COM'CXNUE 

PRINT 91,X(X)/( IJtiD(I,J),UNf,(I,,.T) ,J = 1,MR) 

52 OONTIWIIK 

CAIOi OUTPUT 
STOP 

100 DU 62 T.a1 /#X 

DO 61 0=1, NY 
01(1,0)31.1(1/0) 
r(i,j)3U(i, jj^^'Ud/ J) 

6.1 CONTI tiUE 

62 CONTINUE 

60 CONTINUE 

STOP 
END ■ 



o n n Cl 


'i* 4 *■ 


c 'J.'Hi:y «iiiUfnii-iM.rjr.; CAl.(!iibAXKR kiki,o irri'fcJGiui. 

S 1 1 B K I ) U 'I' X fq (i; X N T K N ( A T N X G ) 

INXI^GKU 

Cin 4 MGiM/(ii/S,p,MX,NR 

C:ilt>M0!'!/().?/X(40) , Ydo) ,H( 10) ,WP,8,l.)JiLX, yw 
CUMMUm/im/FCAOjIO) ,A( 40, 10) ,YBPb(40) ,yBMN{40) 

NYsNH 

SUMaO.O 

Cal. /4./3, 141592654 
Xl = X(fO«XC U+0.5«nE6X 
DU 10 Xai.NX 
X2 = Xl'"lHi;rjX 
YIkUBX'B 

Y2=(y U')-_2,4‘HCn)4!B 
YJaYCf'jVB 

yAaf.YCPJ + '^.’t'HClD^B 

All SATAN (Yl/X2)"ATAM(Yl/Xl)«ATAN(Y3XX2)<-ATAN(y3/Xl) 
A12a«ATAN(Y2/X2)+ATANl y2/Xl)+ATANC Y4/X2)-ATAN(Y4/Xl) 
ACl,t.l=-C*(AIl+AI2) 

SUMs5lll‘i + F(i;,l)*A(I.,l) 

DU 12 J=2#NY 

AHa-A;i2 

Y2 = y2“2,»H( UXtB 

X 4 sY 4 t 2 ,<(H(iJ)«B 

A12a'"ATAN(y2/X2)+ATAN(Y2/Xl)+ATANCY4/X2)-ATAN(Y4/Xl) 
All .iiJ=»C>('CAIl + AI2) 

SlIMsBUM + FCI, Jl^At J. ,J) 

12 CUWTXNUK 

XlaX2 ' ' 

10 cnimwiiK 

A I.A’TrJsfMIM 
RETI.1RW 
END 


THIS Sl.lBRUMTINE BVAMIATES CUEFFICIENXS B'OR WADL INTEGRAL 

subroutiwr: wall 
integer s,p 

CafiMON/Ul/S,P,NX/ NR 
CUMMnN/U2/Gl (40, 10) |W(40,40, 10) 

CfJMMnN/03/X (40 ) , Y C 1 0 ) , H ( 1 0 ) , WP , B , UKDX , Y W 
CaMM0N/U4/F( 40 , 10 ) , A (40 , 1 0) , YBPL(40 ) , YBMN (40 ) 

NYsNR 

CWS"1,/2./3. 141592654 

YlcB»(YCP)-YW) 

yRaB^CYCPl+YW) 



n n o n n n n 


<• 




vyi=Yi*'/i 

'!(X2siy2^Y2 

X).KX(S)-x(U+y.5»nKhx 

XXIaXl^XI. 

AlsAIiOGC CXXl + YYn»CXXl + YX2)) 

FU=A1’AiM(Xl/YU-ATANCXl/y2) 

DU 10 Isi,Mx 

X2aXl«I>Kr.X 

XX2=X2*X2 

a^saIjUgc (XX2+yyj )»cxx 2 tyY 2 ) ) 

B2!sATAW(X2/yn->AXAN(X2/Y2) 

DT = 0.btWP/H»(A2-A.l)-CB2-Bi) 

WC:i,vS,P)aCWiPDl 

Xi=X2 

A1 = A2 

BiaH2 

10 cuNxiMuei 
KKTUHM 
filND 
C 

kid ^ \1» kid kid lA/ Kid pX* kid kL/ kid ^ k.Li ^ ii# tA# di dx «ft# dr dr dr tJr ^ ^ kMM dr dr dU ^ dr d* dr 4W dr df d# df 

THIS SUBROUTINE CALCULATES WALD INTEGRAL 
SUBROUTINE INTWACAIWL) 

Integer s,p 
ccimmon/ui/s,p,nx,nr 

CnHNON/a2/Ul(40, 10) ,w(40,40, 10) 

NyaNR ■ 

5UMr(),0 
DO 10 Ial,NX 

SUMaSUM+irl (I ,Ny)=»'W(I ,S,P) 

0 CUNTIMUPI 
AIWL='SUM 
RETURN 
END 


)(c* )({♦»♦» If 

THIS SUBROUTINE CALCULATES PRESSURE DISTRIBUTION 
SUBROUTINE OUTPUT 

REAL m : ■ I 

, GaMMON/OU/ULLC40f 10),UNL (40,10 ),H,G 
C0MM0N/01/S,P,NX,NR 

COM MON /03 /X ( 40 ) , Y C 1 0 3 , f K 1 0 ) , WP , B , OELX , Y W 
CPC02)aAlf(Cl.+A2’t‘Ci.-Q2))*»EX«l,) , 

92 FORMATdOX^'MAGH NO, a'F7.3,10X, 'PARABOLIC AIRFOIL') 

93 FOBMAT(8X,'X CP LIN CP PRESENT') 

94 FORMAT(bX,F7,3,8Fli;3) 



O C rs cm rs n n a n Cl n n a n a r: n fn n 


m gM 


PHjj4r 
fRlUT 93 
SMsMjkm 
A J,=:‘.^,/(,‘/SM 

f?::x»K/c (}•!,) 

DO 10 laIjNX 

XraXCX) 

DyXRDYCXI) 

UhKCl , +1 }IjIj( K , 1 ) )*»2+DyX#DyX 
aw a ( 1 , + U i\|I:j ( I , i ) ) ♦ ♦ 2 +Dy X ♦DY X 
CPr.sCP(QL) 

CPNaCPCQM) 

PKli'IT 94,XI,CPL,CPN 
10 C 0 f’/'r.TWUB 
RKTURW 
EIWO 


4= 4= If: ^ ^ ^ it: ^ ^ ^ f ^ ^ It it; f )|oK ^ IK * ^ ^ * 4: lie !(! ’(( It! It >1^ ¥ >|( Y 4! 4: if: IK t. !|i ift ^ 

4; 4: 4: 4: i|t if 4: if 4; 4t 4: it: 4: if « 4: 4: 4t :(< it 4: 4: 4; 4; 4: It « itt « « 4( 4! f 4: « it 4; % iit !|: « « « 4t « * it« 4 it « « « !|t « 4( ♦ * It: if « iK « it; « 

THK FDWjflWiNG PACKAGE SHOUIjD BE REPLACED BY APPROPIATE FUNCTIONS 

IF tiie:; body profile is changed 


4c4t4t4t 4: 4<4:4t4(4:4t4(4t 4:414:4: 4:4:44:4:4: 4: 4t4t4t4t4t4:4:4:4<4:4t 4: ittiitttitititiititititiltitVititityitititititicttic 4^4: 4:4: 
THE FULDOWING PACKAGE CALCULATES FUNCTIONS FOR BODY PROFILE 
ITS DERIVATIVE AND LINEAR PBJRTURBATiON VELOCITY FOR A PARABOLIC 
AIRFOIL, 


THIS FUNCTION CALCULATES AIRFOIL PROFILE 

FUNCTION YYtXl 

COMMOM/TY/T 

■■YY«0.,0.:;.. 

lFfXvLT,0,D.CJR,X,GT, 1 . 0 ; return 
YY»2,04:T4X4CI,0"X) 

Return 

SND 


THIS FUMCTION CALCULATES THE DERIVATIVE OF AIRFOIL PROFILE 
FUNCTION DYCX) 

COMMON/TY/T 

OYaO,0 

IFCX.LT.O, 0, OR, X,GT, 1,0) RETURN 
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KfiTURN 

f-lMD 

C 

C 

C THIS rUNCTIUN CALCUr.ATKS I^INKAR PKRl'UKBATIUM VELOCITY 

FUNCTION ULCX^YJ 

COM MOW /U3 / XT ( 40 3 , YT ( 1 0 ) , H ( 1 0 ) , WP , B , DEtiX , YW 
common/ty/t 
PlBi,141592f)54 

lfLs2, O^^T/PI/R^C (0,5~X)+AIjOa( (X*X+Y4‘Y)/( Cl.O 
l-X)i«^2 + Y*Y)3+'4!«0-2,0>KX»CATAN(X/Y)+ArAM((1.0"X3/Y))) 
IF(X.GT,0,.AND»X,1'T,1.) UL=UL+(l,/»»T’t'T>l'(3./PI/PI 
l’t'(2,«2,=*=(X-,S3=»ALUGCX/<U'-X)) 3*»2«1./PI/PI»CAU0G(X/(1,-*X 
2 ) ) 3 >1' il'V.- ( 1 . "4 , =♦< ( X- . 5 ) =t'»2 3 3 ) 

RKTUHrv) 

END 

C 

C 

(J 4! ^ ^ )(t IJ! ;|c »)((!»!))( 4! ^ )(C ♦*:»!;(< !4! !f )((!(( ♦If ^t)|C 
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c , 


;4^ >K f H’ t< ¥ -t ¥ ^ ¥¥ ¥ ¥¥¥ ^ ^ -h -t ¥ -r v ' t: ¥-t ¥ ¥ -r -r ^.. ^ 
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V r Hr’- 
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•C Aiu) / i 

■ '(/.. 5'i'' 

1 1 ; \ 

i-y.' ■ i'i- 



ip< 

f iii'i c riR) i'lU ’■'' 

/ \ 

’V/ 

i<A}- ' 

';i 

} i"' ; ■ 



e'4 
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5 4 i lU'. 

^ .!j 
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\ i'-p;' 
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J * V‘^ i',' 
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i.'itiS’"*] IP ’:itf 7 ';i 

7 y =;(j, , 

( .V , m'(v, O', ; i , , ;s , (f ( ,1,0 i !(i'-'i.'0,%j';' ;■: 

)f > ^■vj0i< I* I'K'J + A + ^ ,0 i; + A. !# 1 h * t A* i. b -A + BH ♦ :< '.J ') ) 
y y-i Ai: ^ " 

feHf- . " 


'n-ilij FUniC'i'lOii CAIjCULiATK^V IVifilcI ^Ari Vii, , t.ii"'; Alj'roib P^OKfLi;:; 

i‘'uwcyi0i4 uycxi' ■ : 

C (71*^ ;AU''^ /'i 1 /T , » Ni.V» B 1 f ij'<i » 1.5 0. / jiOj f B i* , 1 X , : ; . 

py»u,o ■ P P..'',.\ .'//■■■;.•: ■. /■; -'O •■' 

I F cx « t/T', a , C .(j'R.'x, i>T, A vO yiiiulUHll ' > . ■ ■ 

, 0 y so , sip i /sg Ht VX J + b 2 + X * cy 

^ OfisoxW ' ■'.• ■'/■:.■:■■ ^ ■. - 

^'RF’ruRN. ■ o:; ''.''v'-'Vv'. ^"'o:'' :o^:?-v : 


THIS FUNCf lOlt: GALCPli.A’I’gS :^iODiFi:Et?;: LiWEAR FERtURlj ATiUJ^Pf glPOC I 
.■\0t-K fHE'BAXHFOIORV' ''' " 

FUWCTiUR 0B{X) : 

COMHOl{/O3/St40>,PClO)iH(10) 

:cdMyip^i.>'Ty /T dio b 2'/ a 4 ,, ,&6.i a P;,:Y'X ■ 




- i * » I A ) 

U iA ( I ^ U 'i H j ) / lA o t.v'V'C )■' u f { .^ X - v i ^ y :X / i •'•/ i' v'.» 1 > (j 

i'; < M \. 


\ ...; • V •'■ • . ■ •• ■• , ■ 

f U>vL I i H J • C ' jU »C i.U.' ^Vi' ; 1 i ^ Kf ; h MV iv A : ^ r'1 ; i M M WV'lM I F A • ■ sA ' V- MV V ' 

^ A V. \ :Vi! tVKiV (w. Viily. M ^ A Aiij i.j' 

. Mm ' ju ^ / 'V Y /: i’'f V K!.' f ^ V u M t ^-V':’ i* m V' f rvi'i ^ V /•, 

(;i:;^i'.u;:.;f/n 'I /:,. t 4 i,i J » ('■ i, 'I u ; , n U. 0 ) , r , Im-A.;,. 

X i '■ / i . 1. iVi t:V> :i ■: 

• '^'1 C I J \ 

* )*’ I V , i < '!' « c; i j !\ ) 'l i i '4 0 

,. i A»;i: A (X) 

tX , Y ,u *u ) i X . v/M CX I A , f :V;» .i n^' t 

, : ' ■ 

. UQ iij Y , ■ : :, ;,^,' 

■■ Kii=i;:2 + UiXH ■ 

,t-'2«feltt)i.':!i. ' ■ 

S ihVi-'Ri , 4xx^'X tX , )'!l ) I ’I , -l‘j;’*ClX I, ^ V » i'iL* » ^ 

■ : &UM«SUN'+S,lHe ‘ 

10 ■/;. CONTli^Ui:; ", . :■ xx ' 
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